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I. INTRODUCTION 


This thesis is a study of a radiative heating parameterization for 
use in the Fleet Numerical Weather Central (FNWC) prediction system. 
The study has as a primary objective, the analysis of the radiational 
and heat balance of the ocean-atmosphere system utilizing FNWC gridded 
data-fields at constant pressure levels for 16 July 1974. The grid- 
point data were selected along four oceanic meridians, three in the 
Pacific and one in the Atlantic, with most of the gridpoints located 
in the Northern Hemisphere (Fig. 1). 

The specification of amounts of clouds in two designated layers 
has the most influence on the radiative-model dispositions (short- and 
long-wave). The initial specification of the fractional amounts of 
CL(1) and CL(2) were based on large-scale formulations developed by 
Smagorinsky (1960) and used in a similar study based on the data day 
1G January 1974 by Spaeth (1975). 

The governing equations of the radiational transfer of the model 
have been derived by Martin (1972, 1974), who modified both the solar 
and terrestrial radiation transfers to respond to the presence of 
clouds. These clouds were of specified amounts in the two layers, one 
of which is a mid-level cloud and the other a low-level cloud. The 
radiational model in use here has similarities to those in use in UCLA 
and NCAR General Circulation models. 

The physical description of the radiative model may be applied to 


any scale of analysis for which there is adequate resolution of the 
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temperature and moisture data in the vertical. In the horizontal, the 
reliability of the data used here is consistent with that of the FNWC 
analysis to gridpoints, and is typically reported to the nearest tenth 
of a degree with regard to temperature and dewpoint. The radiational 
computations made here are applied to FNWC gridpoints and are designed 
to make a one-hour forward-time a Sagee i Tei to the FNWC primitive 
equation forecast model, with special adaptations to their O-levels. 

The Smagorinsky (1960) cloud-specification of CL(1) and CL(2) when 
used with the model of this study resulted in global albedo estimates 
that were too high especially in the tropics as compared to the satel- 
lite data of Raschke (1973), thus leading to under-estimates of the net 
radiative Suninee at the top of the atmosphere and at the ocean surface. 
This could be a result of inadequately resolved relative humidities, 
which in reality might well be lower than inferred by the FNWC tropical 
analysis. 

The comparative results afforded by the 2/3-CL parameterization 
gave reasonably close agreement with the radiative climatology of 
Peaeedne Ct al (1973), for the NIMBUS III period 16-31 July 1969. 
Oceanic and tropospheric balances were computed for both sets of cloud 
conditions (full-CL and 2/3-CL) for 16 July 1974. It was also neces- 
Sary to augment the purely radiative model by including a turbulent 
boundary-layer model for sensible and latent heat transports over the 
ocean as adapted from Kaitala (1974). 

The radiation package results computed with the 2/3-CL cloud model 
compared very favorably with satellite climatology, though only one 
data day (16 July 1974) was used as representative of mid-July sound- 


ings over the ocean, both for the model and the Raschke satellite 
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climatology. However, the latter climatology was based on the data 
period 16-31 July 1969. 

The major effect noted by the use of the 2/3-CL as contrasted with 
full-CL parameterization was the reduction in global albedo and the 
resulting greater contribution to the surface net heating rates. This 
result is in general agreement with recent observations by Von der Haar 
and Hanson (1969), Von der Haar and Oort (1973) and Raschke et al 
(1973). However, the major discrepancies in the radiative model re- 
sulted from the specification of layered cloud amounts in the tropics 
(in a manner identical to that used poleward of the tropics) where 
the model gave somewhat excessive albedos. It is felt that the geo- 
metry of solar radiation streams impinging upon subgrid sized cumulus 
cells in the tropics should be remodeled so as to divert a greater per- 
centage of solar radiation downward, as compared with the characteris- 
tic cloud-layering and resultant reflective effects which seem 


permissible in middle and high latitudes. 
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II. DATA PREPARATION 


Pee NveEriAL DATA FIELDS 

The temperature and humidity data used in this study were arranged 
in the form of soundings taken along four oceanic meridians (Fig. 1) 
of the Fleet Numerical Weather Central (FNWC) Northern Hemisphere ana- 
lyses on 16 July 1974. Computation of radiational dispositions and 
of other heat budget terms were made at these gridpoints. Oceanic 
locations for these computations were chosen because: 

1. The constant O surfaces (where O = =) of the FNWC primitive 
equation system are close to being constant pressure levels. 

2. The maritime-area heating rate computations are likely to be 
representative of the month of July 1974 as contrasted with computa- 
tions made for 16 July 1974 for a corresponding set of gridpoints 
over land. 

The three meridians (and their respective number of soundings) 
selected over the Pacific Ocean were located at 125W (25 soundings), 
170W (25 soundings) and 145E (17 soundings). The Atlantic Ocean 
meridian was 35W (26 soundings). This method of selecting "sound- 
ings" along the indicated meridians of the FNWC polar stereographic 
map, made it unnecesSary to employ spatial interpolation between 
original data gridpoints along the meridians. Data along line 3 in 
the Pacific was not extended southward of gridpoint (9,55) because 
they fell over land masses (New Guinea and Northern Australia) where 
the surface temperatures and other sounding features were unrepre- 


Snueative OF the Oceanic values. For Similar reason two other 


We, 





(1,63) 


(132) 


(1,1) 


(9,35) 





A= 145°E 


lane 3 





PACIFIC so” 


Eine 2 & (32,32) = A 


now ? 


OCEAN 


Figure l. 


Line 4 
Line 1 
ATLANTIC 
i) > ae CA OC a N 
A=35°W 
A = 125° W 
= (63,1) 


FNWC polar stereographic grid and meridians (lines 1, 2 
3, and 4) selected for study. 


f 
The longitudes A are shown 


for each meridian as well as the extent considered of each 
meridian. 
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gridpoints at (24,24)and (56,6) which also fell just over land masses, 
were replaced with more representative oceanic soundings from the 
erardeoints (23,24)and (56,7), respectively. 

The gridpoint soundings in the form of Table I(a) were taken from 
the original FNWC 63-by-63 surface analysis, as well as the vertical 
distribution of T(P) at ten additional levels ranging up to 100 milli- 
bars. Moisture parameters were provided as the vapor pressure at the 
top of the constant flux level, and as the dewpoint depression at five 
selected levels between the 925 millibar and 400 millibar levels. 

These original soundings were then modeled into a radiative sounding 
(Table I(b)) corresponding in the vertical to FNWC primitive equation 
prediction levels up to 100 mb (See Fig. 2). 

Data from OOOOGMT, 16 July 1974 were used for the Pacific sound- 
ings while 1200GMT, 16 July 1974 data were used for the Atlantic sound- 
ings. These times were used so that the analysis times would correspond 
as closely as possible to local solar noon in each of the indicated 
areas. 

Since the standard instrument level vapor pressure Gel was miss- 
ing from the original FNWC analysis at all gridpoints, Cir 725 approxi- 
mated using the FNWC field of er a computed value of vapor pressure 
at about 20 meters above MSL in the turbulent boundary layer. FNWC 
values of = and e over the ocean had been made available uSing an 
operational planetary boundary layer model detailed by Kaitala (1974). 

To perform radiative calculations, it is necessary to have water- 
vapor and cO., absorber masses and cloud amounts, CL(1) and CL(2) 

Within certain k-level boundaries (Fig. 2). All soundings in this 


study start at sea level with the approximation of surface pressure 


Zi 


Meble ifa). Example of typical FNWC sounding for gridpoint (2,2). 
The humidity parameters between 925, ...,400 mb are dewpoint depres- 
sions, whereas the values at the surface and the top of the constant 
flux layer (9999*) are vapor pressures. 


Pressure (mb) TC) Humidity Parameters 
Surtdce OOO mb 2661 e.. = 23.8 mb 
BS aes bee C 
850 do. Tat 
700 4.8 eZ eit 
500 adds 1355 
400 -23.4 1.3 
300 -38.2 
250 -47.1 
200 = 2)S, 
IES 10 =6/ 2c 
100 -81.0 
9999" te = 23.5 


* Code 9999 indicates data taken from the top of the constant flux 
level (CFL) of the FNWC initial data program. 
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Table I(b). 


and CO 


Smagorinsky 


Example of the corresponding radiative sounding with mix- 
ing ratio listed at odd k-levels (Fig. 2). 


Additionally, water-vapor 


absorber masses, and cloud amounts CL(1) and CL(2) (after 
(1960)) are also listed, as these parameters have been 


modeled in the radiative theory presented in this study. 
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Figure 2. Five-layer radiative sounding used in this study. Levels 

are identified by their values on the k-scale, while layers are identi- 
fied by their level boundary indices in parentheses, e.g. (8,10). 
Pressure-scaled water vapor and CO, mass increments AU and AC, respec- 
tively are integrated relative to the surface and the resulting U and C 
are carried at even levels. The temperature T is retained at all levels. 
Amounts of clouds CL(1) and CL(2) in the layers shown have been para- 
meterized for consideration of their radiative effects. 
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1? = 1000 mb. Therefore, the k-levels are associated approximately with 
the FNWC levels PL =O Oe Omer oUOn, eee, 2UUn, U0”, 0.0 ma respec— 


tive EO a = 10; Oso, e@ey ems O20. 


B. INTERPOLATIVE PROCESSING TO K-LEVELS IN RADIATIVE SOUNDINGS 
1. Temperature Profiles 
The gridpoint temperatures were listed at each mandatory level 
of Table I(a) between 1000., ..., 100 mb. The temperature was assumed to 


be isothermal from 100 mb to p=0.0 mb. The temperature T was set equal 


10 
to the FNWC listed sea-surface temperature. The radiative sounding tem- 
peratures for the remaining k-levels were obtained from either their 
corresponding listed temperature value or by a three-point Lagrangian 
interpolation scheme (Eq. 2-1, (Spaeth, 1975)), to level k when the 
listed temperature-profile does not include the value Ty 
2. Moisture Profiles 

The moisture from the original FNWC soundings (Table I(a)) were 

converted into mixing ratios at each of the original sounding levels. 


The near-surface vapor pressure (taken as ) was used to calculate the 


mixing ratio at k=10, as described by Spaeth (1975): 


qo = (6b Sy (e/1000) (2-1) 


The remaining mixing ratios q were calculated from the original FNWC 
sounding (Table I(a)), using a dewpoint depression formula (Spaeth, 1975, 
(Eq. 2-3)). The computed gq-values were subsequently interpolated to k- 
levels using the previously mentioned three-point Lagrangian procedure 
applied to successive g-values in the original sounding. Resulting 
values of the interpolated g-values are shown for the case of the radia- 


Pivessounaing at Gridpoint (2,2) (Table I(b)). 
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Due to the fact that most radiosonde humidity data for p<300 mb 
is either unreliable or not available, a power-law extrapolation for- 


mula (Spaeth, 19/75) 


: r 
W/o = (P| /500) (2-2) 


was used to obtain q-values at k=3,2,1. 

This extrapolative procedure for obtaining q-values at high 
levels was first suggested by Smith (1966), utilizing the least squares 
estimation of A according to 

6 


i=1 Yi : 


A= —+——— = R (2-3) 





-M 


with a= log (q./q.) and x = leg (D/P .) The Smith method was then 
further tested by Spaeth (1975) for the vertical scale of the FNWC 
initial data soundings. The resulting correlation coefficient R 
was found here to lie in the range .95 to .99. Such high correlations 
indicate that the procedure described by Spaeth (1975) for determining 
upper atmosphere q-values has realistic estimation value for the verti- 
cal scale involved in the radiative sounding. 
3. Pressure-Scaled Absorber Masses 

The pressure-scaled water vapor absorber mass AU ina layer 
(Fig. 2) were calculated for the five odd k-levels by using the computed 
mixing ratio values (Eq. 2-7, Spaeth, (1975)). The equation for com- 
puting the integrated pressure-scaled water vapor mass is given by the 


algorithm described by Eq. 2-8 of Spaeth. 
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A similar routine was used in the computation of the carbon- 
dioxide scaled absorber mass (Spaeth, 1975, Eqs. 2~10,11). To be noted 
mere 1S the fact that cO., absorber masses have been stated in terms of 


N.T.P. “pressure-scaled"cms, since the CO, absorption coefficients are 


Z 


generally stated in terms of the reciprocal of this CO. mass unit. 


eee chOUD PARAMETERIZATION 

The relative humidities and thus the saturation vapor pressure at 
levels k=5 and k=9 are used in the calculations of the fractional cloud 
cover CL, and CL., in layers (4,6) and (8,9) respectively, as depicted 
in Fig. 2. From the computed relative humidity values at k=5 and k=9, 


the fractional cloud amounts for each layer were parameterized using 


the equation (after Smagorinsky, (1960)): 


CL (1) Ze Oeeert (5) = 10. 7 (2-4 (a)) 


Bro OUR 9) 2) peas O (2—4 ())) 


CN 2), 


Smagorinsky's parameterization of CL(1) and CL(2) permits cloudiness 
fractions of 1.30 and 1.333 respectively with RH=1.0. Values of CL 
greater than 1.0 were considered initially by Smagorinsky (1960) to 
suggest the amount of supersaturation which accompanies precipitation. 

All radiative models in use at the present time limit CL in accord- 
amnee with O < CL < 1.0. In the present radiative package there is little 
opportunity to deduce precipitation and/or supersaturation amounts. 


Hence, the formulations of CL. and CL 


1 > of Eqs. 2-4 were limited toa 


Hasamum OL L.O in each layer in the “full-CL" case. In addition, the 


27 





Smagorinsky computations for CL, and CL. were reduced to 


O) 
a 
Tt 


O72 (CL, ) (2-5 (a) ) 


O 
tt 
{| 


2/3 (CL) (2-5 ()) 


regardless of whether the original CL values exceeded 1.0 or not, so 
that the revised cloud-cover amounts could be considered for purposes 
of estimating shifts in the global radiative balance. Computations 
made with cloud parameterization of Eq. 2-5, are termedthe 2/3-CL case. 


gae 1/3 reduction inherent in CL.', CL 


1 a relative to Smagorinsky's 


CL, CL., were initiated in this study in an attempt to tune cloud 
amounts to give radiation results in closer agreement with recent satel- 
lite climatology of Raschke et al, (1973). This satellite climatology 
suggests the use of smaller cloud amounts than that specified by 
Smagorinsky. Fractional cloud amounts were considered to be functions 
of large scale effects only. Therefore, small scale convective acti- 
vity, seasonal conditions (except for the soundings) and latitudinal 
effects were not considered in specifying the reduction factor of 


Eq. 2-5. The reduced cloud-parameterization was introduced here for 


estimating large scale radiational effects only. 


D. CLOUD-AREA COVERAGES 
Since Eqs. 2-4a,b or 2-5a,b gave the fractional-coverage of the 
gridpoint area by the appropriate cloud-type, the gridpoint area may 


be thought of as broken into random fractional segments of size 


Wwi0,0). = (1-CL, )* (1-CL,) (2~6 (a) ) 
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wherein there iS a combination of clear-over-clear segments in the 
layers. Similarly, the gridpoint area has the fractional area of cloud 


coverage 


W(1,1) = CL oC (2-6 (b) ) 


of upper-cloud amount cL, overlying lower-cloud amount CL... Likewise 
the combinations of cloudy-over-clear and clear-over-cloudy, by layers, 


May be visualized as occurring with the weights 


WwW Ge0) 


I 


CL, = (1-CL.) C2-6(e)}7 


and 


W(0O,1) 


(1-CL, ) 2 CL (2-6 (d)) 


Z 


Regardless of whether the full-cloud amounts CL,, CL, of Smagorin- 


1 Z 


Sky (Eq. 2-4), or the 2/3-CL amounts of Eq. 2-5 were utilized, it was 
useful to carry, for radiation computations of each sounding, the 
relative weights or fractions of the gridpoint area exposed to the 
specified cloud-layer combinations. Henceforth, the symbols are denoted 
Pea, 0), W(1,1), W(1,0), W(0,1) and as given by Eqs. 2-6(a,b,c,qd) 
Suggest overcast (1) or clear (0) cloud amounts in the indicated layers 
tetg. 2), the first index referring to the level of CL, and the second 
Perthat of CL.- The usefulness of this computational device will be 
exemplified in Sections III and IV, where the procedures for the terres- 


trial and solar radiational computations are discussed and the results 


are summarized over the full set of soundings. 
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i.) LE RRESTRIAL RADIATION 


A, THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas were developed by Sasamori (1968) for flux emis- 
Sivities in the atmosphere as associated with computations required for 
the heat balance requirement of the NCAR General Circulation Model. 
Sasamori derived the empirical emissivity formulas for water vapor and 
cO., by comparison with the theoretical values built into the Yamamoto 
Radiation Chart (1952). This chart has proved to be quite accurate for 
computational checking of the SaSamori emiSSivities and was used by 
Spaeth (1975) and Warner (1974) as a systematic guide for integration 
of the radiative transfer formula developed by Martin (1972, 1974), who 
adapted the Sasamori formulas to the particular layers of interest in 
gridpoint computations of the FNWC primitive equation model (Fig. 2). 


The essential long-wave flux formulas reguired for use in the FNWC 


heating package are the following: 


Foo = net IR flux at earth, k=10 
* 

Po = net IR flux at level k=6 
* 

F —wiCt eerie at. Level k=Z2 


F610 = net IR-flux divergence in the layer (6,10) 
F26 = net IR-flux divergence in the layer (2,6) 
To get the IR-flux divergences in the layers (6,10) and (2,6), the 


* * x * 
differences F - F and F -~ F must be computed. The detailed 


6 10 2 6 
scheme for making such computations with various combinations of cloud 


cover CL(1) and CL(2) is similar to that developed by Martin (1974) and 


as reproduced in detail by Spaeth (1975). 
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In order to make IR net-flux calculations along the path of integra- 
tion, there must be a physically sound representation of the emissivity 
 ) as a function of both water-vapor and co, absorber masses in layers 
along the sounding. For a further discussion of the emissivity formulas 


used in the quadrature scheme, refer to Spaeth's Appendix A (1975). 


pee NET FLUX FORMULATION 
1. At Level k=10 

The radiative sounding as depicted in Table I(b), was computed 
as the combination of parameters U(k,10), C(k,10) and Ty for each re- 
maared level, k=10,8,...,2,1,0. Cloud parameters CL(1) and CL(2) are 
also listed at each gridpoint and in general are non-zero. The grid 
area was then considered to be composed of areal fractions (weights) 
@efined in Eqs. 2~-6(a,...,d) and denoted by the symbols W(0,0), W(1,1), 
reol,O), W(0,1). 


The composite net flux F (CLs CL.) at level k=10 at each 


LO 


gridpoint is then constructed by using the appropriate weight factor to 


* 
multiply the corresponding reference net flux F computations ,defined 


10 


* * * 
for the special cloud-cover cases>F (O07 sr (le One oF (Giles 


10 10 10 


* 
Foo (1,1) .As a result, the composite net flux may be written 


x * * x 
1S (CL, CL.) = Oe (0,0)+W(1,0)F, (1,0)+W(0,1)F) (0,1) 


10 0 0 


* 
(1,1) (3-1) 


0 


+W(1, 
(1 LF, 5 


Spaeth (1975) has listed these reference net flux formulations in his 


Sul 





fos. 3-6,7,8. Using the definitions of W(0,0), W(1,0), W(0,1), W(1,1) 


* 
19 (CL 


10 prChy) can be shown to assume the form 


* 
F (CL, ,CL.) = [1-CL(2)}1(B, 5B 


10 )~-5{€ (8,10) (B, .-B,) 


6 
+ (€ (8,10)+€ (6,10)) (BL~B.)]} 
+ (1-CL(2)) (1-CL(1)){B_-.5[(e€ (6,10) 
6 wc 
+¢€ (4,10)) (BL-B,)+(e€ (4,10) (3-2) 
wc 6 4 wc 
+¢€ (2,10)) (B.-B_.)+(€ (2,10) 
WC Ae wc 


4 €_ (110)) (BB, ) +E (0,10) ,T,) *B, J} 


+ 


CL(2){(B, .~B,) [1-.5€,,_(9,10)]} ; 


Here 


B, = 1.170403 x 10’ in (geo) 


is the Stefan-Boltzmann blackbody flux in ly AE a at temperature Ty 

Further, EO Ur Cys 10) is the combined water-vapor and cO., 
emissivity along the path from level 10 to level k. This emissivity is 
considered by Sasamori to be temperature independent for T > 210K, 
whereas ey represents the temperature dependent emissivity applicable 
for 7 < 210K (see pp. 136-137, Spaeth, 1975). 

The reference net fluxes Pro. of Eq. 3-1 are associated with (1) 
Clear skies in both layers, (2) overcast in the upper layer only, (3) 


Overcast in the lower layer only and (4) overcast in both layers, 


respectively. 
* 
6—- 


* 
The formula for Po (CL, ,CL,) has been developed by Spaeth (Eq. 


Pee Net Flux F 


3-10, 1975) in a manner analogous to the Gerivation of the weighted 





F ™ fhe result is reproduced for completeness after (pp. 48, Spaeth, 


FPF. = [1-CL(1)]{B,-.5[€ (6,8) (B,-B.)+€, (4,6) (B.-B,) 
+ (6 (4,6)+€ (2,6)) (B,-B,)+(€ (2,6)+€ (1,6) (B,-B)) 
+ € (0,6) ,T,) *B,]}+(1-CL(1)) (1-CL(2)) {(B, (Bg) [1. 
-.5(€ (6,8)+€  (6,10))]}+CL(1){(B.-B.)* 

wc wc B36 
ils eee) tein) (11, (2))) 1 (Bb) --B.) * 

WC LOS 


Per Sioemioscite (Gc 1oV01) . (3-4) 
WC WC 


x 
3. Net Flux F,_ 
The net flux at level k=2 can be calculated in a similar fashion 
* * « ° « 
EO Foo and Fe , with the formulation of the individual reference net 


* * 
fluxes P., as previously described. The formula for F. (CL, ,CL,) is 


* * 
developed analogous to that for F (CL) ,CL.) and Pe. (CL, ,CL.)- 


10 


* 
[= = er Uaioe 7k) SBD) +e (2,4) 
= Eo (216)) (BE-B,)+(€ (2,6) +E, (2,8))* 
(Bo-B)+E (1,2) (BL-B) +e, (0,2), TOME ee 
+ (1-CL(1)) (1-CL(2))*{(B..-B.)[1-.5(e€ (2,8) 

KO Es. wc 
re €_ (2710)) ]}4+cL (1) {B,-.5le (2,4) (B,-B,) 


+ € (1,2) (B,-B)) +e, ((0,2), oc) eae (3-5) 


A new parameter, the total outgoing long-wave radiation to space 


* 
(FF2) can be readily defined from the expression for F by setting to 


2 
zero the terms of Eq. 3-5 representing the downward flux through level 


k=2. These terms are just those denoted by the symbols See One) a) 28, 


a)5 





and Eg (112) *(BL-B)- If the long-wave radiative model including its 
cloud parameterization is realistic, the modeled outgoing flux to space 
(FF2) should compare reasonably with that observed by NIMBUS III satel- 
lite during the same period, after Raschke et al (1973). The most 
meatly comparable satellite climatological period was 16-31 July 1969. 
These satellite observations were obtained from the NIMBUS III satellite 
atlas of Raschke et al (1973) for the July dates indicated and used for 
Bemparison with the model values of FF2, point for point (at 5° latitude 
intervals) along the same four meridians of Fig. 1. Also the latitu- 
dinally-distributed mean satellite fluxes were compared with the corres- 
ponding averaged values of FF2. The results showing the comparisons of 
the latitudinallyzaveraged distribution of the model computations of 

FF2 and the corresponding satellite observations are listed by latitude 


im Table V. 


C. APPLICATIONS TO HEAT BALANCE COMPUTATIONS 
1. At Upper and Lower Boundaries 
In order to compute the heat balance at the top of the earth- 


* 
tropospheric system the composite F and the total insolation absorbed 


2 
below level k=2 at each gridpoint are required. Moreover, a radiative 


balance is computed at the surface for each gridpoint using the earth's 


x 


absorbed insolation (Section IV) and the composite net flux loss Foo 


In order to consider the possibility of a heat balance at the 
earth's surface, rather than merely a radiative balance, it is necessary 
to add a term representing combined sensible and latent heat transfer 


meross the alr-sea interface (Section V). This latter transfer rate at 


each gridpoint was adapted from the FNWC primitive equation model, after 
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the discussion of Kaitala (1974) and based upon the FNWC data of 16 
July 1974, 
2. Intermediate Levels 
Consideration of the radiation balance in the atmosphere re- 
quires computation of the long-wave cooling effects caused by the 


flux divergences F26 and F610. F26 and F610 are defined as 


F26 = F =r (3-6) 


F610 (3-7) 


Il 
Ny 
| 
ry 


where the symbols (2,6) and (6,10) indicate the layer boundaries in- 
volved in the balance consideration for the indicated layer. 


In order to compute F26 and F610 at each gridpoint, complete 


* * 


1S 


* 
10 ' Fe ' Fos associated with each set of 


listings of IR net fluxes F 
mererence cloud amounts (0,0), (1,0), (0,1), (1,1) have been computed 


* 
at each gridpoint together with the weighted set of fluxes Foo (CL, 


* * 


oy, ©, (CL,, CL), F, 


6 1 ECL 


1’ CL.)- These cloud-weighted values have 
been constructed using the weighting scheme of Eqs. 2-6, 3-1, 3-4 and 
3-5. A sample gridpoint printout of the net fluxes and flux divergences 
has been included in Table II. 

The determination of the heat balance by layers has been de- 


ferred to Section VI, where the results are displayed in cross- 


sectional form for the 2/3-CL case, 


Pee OLATISTICAL RESULTS AND COMPARISONS 


* 
1. Net Flux Foo CO. 0) 


* 
A statistical test of the Foo (0,0) numerical results was con- 


ducted by utilizing a linear regression program for relating the 
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* 
(0,0) against the simultaneous values of X. = B and 


predictand F 1 10 


19, 


x, = male as predictors, uSing all 93 gridpoint data as samples. The 
regression program was based on the BMDO2R in the Biomedical set of 


programs (Dixon (1973)). The result was cast in the form of the well- 


known Brunt net flux equation (1932): 


hr} 
i 


- B,  (atbve) (B28) 


where 


= * = 
Blo St Tio (3-9) 


is the Stefan-Boltzmann blackbody flux at the surface temperature Tio 
and e is the surface vapor pressure in mb. The multiple correlation 
Secerficient RUE ere) was higher than 0.99 in the 16 July case, 
as was also true for the data periods of Warner (16 October 1973) and 


Meyers (16 April 1974), (Table III). 


* 
Since Foo may be defined as 


F = 5 - F (3-10) 


fem o-lLO can then be solved for Fa: 


Pa = 8 Gl bve] (se119 


for the clear-sky case. It must be recalled that the computation of 


Fs included both the effects of water vapor and of CO Since there is 


De 
no direct relationship of Ye upon the quaSi-constant absorber mass of 


Con, the radiative effect of the cO., in Eq. 3-11 must be included in 


the term involving (l-a) Bio: 
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A cross-Seasonal comparison of the coefficients l-a and -b of 
Eq. 3-11 is presented in Table III, comparing the clear-sky results of 


Fs for 16 July 1975 with those of Meyers and Warner at climatological in- 


tervals of three months earlier and later, respectively. 
The high multiple correlation coefficients indicate that an 
oceanic version of the Brunt downward flux equation is valid for each 


tested season. Fs depends most strongly upon Bio while only slightly on 


B ave: This is an apparent result of the nearly constant relative 
humidity over the large oceanic regimes. Moreover, the weak dependence 


of F_, upon B Ye is most evident at the data-time 16 July 1974, when 


d 10 


the coefficient -b is so small as to be statistically insignificant. 
This result is due primarily to the relatively weak gradients of both 
Ye and in the radiative soundings aloft corresponding to 16 July 1974, 
as compared with the adjacent seasons. 


* 
Zee Modification of Pio for Cloudiness CLis CL. 


A second statistical test was performed based on a possible rela- 


* 


10° The 


tionship between the cloudy-sky and the clear-sky cases of F 


* * 
ratio Y = F (CL, ,CL,)/F, (0,0) was used as the predictand in this 


10 O 


case. The predictor CL was the total opaque cloud cover at each grid- 


Pent, defined as (after Quinn, 1971,and Spaeth, 1975): 
Gig = Chit ye 406 (2) (12=CL (1) ) (3-12) 


* 
This regression was run utilizing F (CL,, CL.) in the numerator, com- 


10 


puted using the amounts of CL(1) and CL(2) as calculated from Smagorin- 


* 
sky's formulation and then utilizing F (CL, ', cL,') computed from the 


ae) 
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2/3-CL parameterization. The general form of the regression formula 


required is shown in Eq. 3-13 below, with CL given in Eq. 3-12 


FP) (CL) == Fy (0,0) (1-a(CL) ) (3-13) 
A cross-seasonal comparison of the coefficient d for the 2/3-CL and 
full-CL cases using the results of this study and those of Meyers and 
Warner, and of the resulting multiple correlation coefficients are pre- 
sented in Table IV. There seems to be little clear-cut variation in 
the statistics of Table IV between the results of adjacent seasons, 
and little statistical preference between the full-CL and the 2/3-CL 
cases. The high correlation coefficients of Table IV, for all seasons 
considered, indicate the general capability of the IR radiative model 
employed in this study to account for the effective net radiation at 
the surface. 
k 
The mean values of F over the 93 gridpoint soundings for 


10 
16 July 1974, as applicable in Tables III and IV were as follows: 


PVli—CheStatistics 2/3s—CiaSstatistics 

= Se eck a Soe 

tee (Ci) = .0774 ly min Ff.  *(CL )=. .0966 ly min 
10 10 

SSS wel a toed 

Poe O,0)— 1449 ly min Po 2 (0,0) = 21449 1y min 
10 tO 

CL = .5890 CL'= .4306 


The superior bar symbol denotes the sample mean. These results indi- 
cate that the surface net flux is decreased by the ratio .4658 with 
CL = PoovO, and by the ratic ~3333 with CL' = .4306. These results 
indicate that the downward flux F (CL) ESuOvoupercent Larger at the 
Subtace than P ,0,0) when full-CL amounts are used in the cloud para- 


meterization model. In the 2/3-CL parameterization the mean downward 
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oer ilh 1974 Vo July 1974 6 Occover, 1972 


ea .6557 .8129 6834 

~b .0209 -.0094 .0161 
re .4527 .5471 .4716 
Ry .9969 .9950 .9973 


Table III. Cross-seasonal comparisons of the coefficients l-a and -b 
from Eq. 3-11, as well as mean values from the Brunt equation and the 
resulting multiple correlation coefficients have been included. 


16 April 1974 16 July 1974 16 October 1973 
2/3-CL full-CL 2/3-CL full-CL 2/3-CL full-CL 
d = 8 375 -./4 soot oll wed TLE =a 
R -I9 16 ~9945 s9942 72920 9974 soo72 


Table IV. Cross-seasonal comparison of the coefficient d from 
Eq. 3-13 and the resulting multiple correlation coefficients for both 
the 2/3-CL and the full-CL cases. 
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miux at the surface PfcL’) is increased by only 33.3 percent when 


compared to P 60,0). 


E. COMPARISONS OF FF2 WITH SATELLITE CLIMATOLOGY 

Comparison was made of computed-model values of FF2 with satellite 
measurements of total long-wave flux to space (after Raschke et al, 
M73) for the NIMBUS III period of 16-31 July 1969. FF2 is used here 
instead of F, Since the downward flux at k=2 is not measured by the 
satellite. The computation of FF2 was previously discussed. FF2 for 
16 July 1974 was computed for both full-CL and 2/3-CL cases at all 93 
gridpoints. These FF2 values were then interpolated to whole multiples 
of 5° latitude along each of the four meridians and the results were 
then averaged zonally to get a mean latitudinal distribution valid for 
each 5° latitude interval considered in this study. 

Table V shows both of the latitudinally-~distributed FF2 cloud- 
model values compared with those extracted from Raschke et al (1973). 
Raschke's results were obtained by interpolation from charts using the 
same oceanic meridians as those used in this study. Latitudinally dis- 
tributed total opaque cloud cover CL by Eq. 3-12 are also listed for 
both full-CL and 2/3-CL cases at each latitude. In the bottom line 
of each column in Table V is listed the cosine-weighted mean of each 
set of the column values for the listed latitude range. It should be 
noted that the Raschke results do not correspond to specifically 
known CL values. Table V suggests that the model results for FF2, 
both full-CL and 2/3-CL, are reasonably close to the corresponding 
values derived from Raschke. A close examination of this table shows 


that in the regions 20S to 5N and 25N to 45N both sets of model results 
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Table V. Comparison of zonally-averaged longwave flux to space, FF2, 

aS found by this study for both full-CL and 2/3~-CL cloud cases for 

16 July 1974, and by Raschke et al (1973) based upon NIMBUS III measure- 
ments. Also included are composite cloud amount fractions used in the 
evo cloud models of this study. Flux values in ly min71, 





are slightly less than given by Raschke, while in the region 10N to 
20N both sets of values are ene higher than Raschke. From 20S 

to 45N the agreement between FF2 and F2(RAS) is quite close. From 

50N and higher the correspondence to F2(RAS) is not as good and perhaps 
this is due to the high amounts of cloud cover experienced on this 

Gey an this region. This could be particularly true since the model 
calculations are for a Single day, whereas Raschke's calculations 

are based upon 16-day averaged values. Additional limitations of the 
comparisons made here are obvious, when it is recalled that between 
latitudes 20S-5S and between 60N-65N there are fewer than four meri- 
dional lines available for computing the listed zonal values in 

Table V. For all other latitudinal~average values, four meridional 
lines were used in the averaging. While the results from Table V are 
not used to make a final choice of cloud-cover parameterization for 
use in the radiation model, there is evidence of slightly closer agree- 


ment of FF2 with F2(RAS) using the 2/3-CL parameterization. 


F. COMPARISONS OF CROSS-SEASONAL RESULTS OF Plo 


* 
Table VI depicts the latitudinally distributed values of Fao (CL') 


and of CL‘ at the earth's surface obtained using the 2/3-CL parameteri- 


* 


zation. A seasonal comparison of the model computed Foo is made with 


data obtained from Meyers (1975) for the 16 April 1974 case and 


fmm Warner (1974) for 16 October 1973. It is clearly shown that 


* 


Foo (CL') is a decreasing function of cloud cover. The most noteworthy 


c 
aspect of Table VI is the belt of high values of F fom to —Jullyy19 74 


10 
memeene latitudinal band 20S to 10S. If the associated CL values are 


examined, it is noted that the composite cloud values were very low 
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Table VI. Cross-seasonal comparison of the zonally-averaged net 
mluxes at the surface F i and zonally-averaged CL‘ based on 2/3- 
CL parameterization for the three mid-seasonal dates 16 April 1974, 
moewuly 1974, and 16 October 1973. Flux values are in ly Manne 
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and this seems to be consistent with the high values for F Further- 


* 
Orr 
more, the low CL values would appear to be consistent with the northward 
migration of the ITCZ during mid-summer. There is a tendency in each 

* 
season for a relative maximum value of Pio to be located in the sub- 
tropics (latitudes 15-251). There is further evidence of a high latitude 


* 
(55-65N) minimum F associated with a concentration of maximum cloud 


10 
cover CL'. The analysis does seem to be consistent across seasons. Simi- 
lar conclusions were reached for the full-CL model, although they were 
mee iunciuded in Table VI. In general, it is felt that the results of 
the ae analysis yield no conclusive results in regard to the appropri- 


ateness of the choice of full-CL or the 2/3-CL parameterization models. 
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iV. SOLAR RADIATION 


Ae PARTITION OF SOLAR INSOLATION 

The solar constant assumed in this study at level k=O (top of atmos- 
phere) was 2.00 ly Pere (Joseph, 1971). This flux was depleted by 43% 
to account for the attenuation caused by oxygen and oZone above the tropo- 
pause. This left the value S = 1.92 ly age at level k=2 to be used 
in this study as the effective solar constant. 

Equation (4-1) was then utilized to compute the effective solar in- 


solation at the tropopause (k=2) as follows: 


BUD = & (= Cos z (A213) 


m 


II 


where S effective solar constant at k=2 
Cos z = cosine of the zenith angle for the Julian date used 


x/x = ratio of the actual earth-sun distance to the mean earth-sun 
distance for the Julian date used in this study. 


The Smithsonian Meteorological Tables (List, 1958) gives the ratio r/X_ 
and the solar declination 6 for 16 July 1974, OOOOGMT, as these values 


are employed in Eq. 4-1. These values are, respectively: 


II 


1.01644 


Il 


Yr 
ze 
m 
6 


21.5N latitude 
6 is used in evaluating the cosine of the solar zenith angle, given by 


Cos z = Sin $ Sin 6 + Cos ¢ Cos 6 Cos h (4-2) 
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where @g@ is the latitude and h is the hour angle of the sun relative 
to the meridional data lines. For example, Fig. 1 makes it clear that 
for 

line l h = 55° 

line 2 lin = 20:7 


55. 


lines 3,4 h 
at the times of the synoptic charts (OOOOGMT and 1200GMT, as applicable). 
Sin $ waS computed according to the standard polar stereographic projec- 


tion formula applicable to the FNWC base chart given by 


rae ~ £(1-32)7 + (3-32)73 


Sin 6 = (4=3) 


a ste [ (I-32) 7 che a3 a] 


OE ' 
where re = 973.752. Thus Sin ¢ assumes the following functional form 
in terms of the FNWC grid-coordinate I (Fig. 1): 
ines 1, 3, 4 


2 
Sin ¢ = 2232752 =2(32-1) VAR 


2 
SIG ela a2 (52-2) 


Erne 2 


2 
Sin @ = 223+752 = (32-1) aay 


2 
Oe Sows ee (3 2=1) 
or, conversely I is given in terms of ¢ by 


mines 1, 3, 4 


z = _cos o> 2 
I = 32 - 22.065[>——Zo aI (4~4(c )) 
Line 2 1, 9D ee oe ] (4-4 (a)) 


Wades 1 


ieee co Or Lines 1, 2 
[eos oo OL bane 3 


ij OS. oo. FOr Line.4 


4/7 





For lines 1, 2, and 3 the gridpoint soundings correspond to OOOOGMT 16 
July when solar noon occurs at the 180th meridian. Line 4 gridpoint 
soundings correspond to 1200GMT, 16 July, when solar noon was at the 
Greenwich meridian. 

A very simple partition of solar insolation was utilized in this 
study after Joseph (1971). It consisted of dividing the insolation F(2) 
into two parts at level k=2. One part was considered to include all 
wavelengths A > .9 um where absorption by water vapor and carbon 
dioxide bands are the most prevalent attenuation processes in clear air. 
This part of the solar spectrum was termed F(A) energy and was considered 
subject to water-vapor absorption but not to Rayleigh scattering. For 
those wavelengths A fe] ico semP rien GEathe selar insolatien energy 
by water vapor was considered negligible. This part of the solar insola- 
tion was denoted F(S) suggestive of the fact that it was subject only to 
Rayleigh scattering attenuation in clear air. The two partitions are 


formulated after Joseph (1971) as follows: 


F(A) ~349 F(2) (4~5(a)) 


I 


tS) FOS ParuiZ) (4=5'(5) ) 


II 


In this study, the introduction of two cloud decks produced cloud-reflec- 
tivity effects upon both the F(A) and F(S) solar energy insolations. 

However, in the clear areas around any gridpoint the absorption-attenua- 
tion only applies to the F(A) insolation, while only Rayleigh scattering- 


attenuation applies to the F(S) insolation. 


Pee DISPOSITION OF F(S) INSOLATION 
In the disposition of the F(S) insolation, Joseph (1971) determined 


that Rayleigh scattering reflectance to space by clear skies (after 
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Coulson, 1959) could be effectively approximated by least squares in 


the following form 


a(R) = .085 + .25074 [log (=— Sec z)) (ae) 


Oo 


where Eee 1013.25 mb. In Equation 4-6, T/P —iermnewiew Of the fact 


that the mean sea level pressure TT is close to 1000 mb. Also 
-1 
Sec z = (Cos 2) 


weasel COS Z given by Equation 4-2, 
The surface albedo a(G) is another reflective parameter utilized 
in this study. Over oceanic areas the following formula for Q(G) 


after Gates et al (1971), was utilized: 


a(G) = max {.06, .06 + .54 (.7 - Cos z)}. (4-7) 


in Clear Sky Case 
In the clear sky (0,0) case the F(S) insolation was subjected 
to Rayleigh scattering-reflectance a(R) by the atmosphere, as well 
as to reflectance a(R) by the earth's surface. Considering the likeli- 
hood of a succession of multiple reflections between earth and atmosphere, 
the F(S) insolation actually penetrating the earth's surface after 


scattering is given by 


IS10(0,0) = F(S) [1-a(R)]} [140 (R)a(G)+.... (a(R)A(G))™ 
+..0)* (1-0 (G)) (4-8 (a)) 


Baeat 1S, by 


Oe O) meets ot) | [i-O(G) ] / (1-a (Ry a(G) J (4-8 (b) ) 
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Pe Loud y—-ony Cases 
In the three cases in which clouds were present, F(S) insola- 


tion absorbed by the ground at each gridpoint was (Arakawa, 1972): 


TS10(1,1) = F(S) (1-R(1)) (1-R(2)) (1-a(G) ) 
*{1-[(R(1)R(2) + R(2)a(G) + R(1)a(G) 


+ 2R(1)R(2)a(G)}} 2 (4-9) 


memindicated by the notation (1,1), denoting CL(1) = CL(2) = 1.0, 
Eq. 4-9 is the formula used in calculating F(S) insolation absorbed by 
the ground in the case where overcast clouds are present at both levels 
of Fig. 2. Also in Eq. 4-9, constant cloud-reflectance values were 
chosen, namely R(1) = .54 for the mid-level clouds between k=4 and 6, 
and R(2) = .66 for the low-level clouds between k=8 and 9. Both cloud- 
reflectance values are as suggested by C. D. Rodgers (1967). 

For all the other cloud cases, the following changes were applied 
moera.e 4-9, In the (1,0) case (CL(1) = 1.0, CL(2) = 0.0), the desired 
result is obtained by setting R(2) = 0 in (4-9), from which it follows 


=nat 
oO (dO) P= eb (5) Cl-R(1) ) (faa (G)) 7 [1-R(1) (6G) J (4-10) 


In the case (0,1), one sets R(1) = 0.0 in (4-9) so that (4-9) is simpli- 


fied to 
PSOKO nL) 9= bis) (1-R(2)) (1-@(G) ) / (1L-R(2)a(G) ] (4-11) 


Note that with a cloud overcast present, the Rayleigh clear-sky scatter- 
ing a(R) does not appear in Eqs. 4-9, 4-10 or 4-11, but is included 


in the cloud reflectances (Rogers, 1967). 
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ce Composite E(S) insolation 


Equations (4-8), (4-9), (4-10), and (4-11) were utilized in the 
computation of the cloud-weighted F(S) insolation penetrating the 
earth's surface considering the areal-weights of the cloud combinations 
Memoced by (0,0), (1,1), (1,0) and (0,1) about a gridpoint.. The resultant 
F(S) insolation penetrating the earth's surface denoted by IS10 is there- 


fore expressible as 


noone Cpe Ch (2) = ESTO(O.10) WwiO,0) 


ob 


TS TOs lea a2. 1) 


+ 


PeLO<d, Oj wl, 0) 


+ 


IS10(0,1) W(0,1) (4-12) 


Here the weighting factors W(0,0), W(1,1), W(1,0) and W(0,1) are as com- 
puted in Eq. 2-G(a,...,d) respectively. Note that Table VII lists the 
radiative computations as they apply to the disposition of F(S) insola- 
tion. The four individual computations of IS10 as applicable to the 
possible overcast-clear layer cases are made under that heading while 
the composite value of IS10 appears in the final line of the same column 


Of Table VII. The difference 
REFS = F(S) - IS1O (4-13) 


in each case represents F(S)-insolation reflection to space while 


. IS10 _ 
STRAN = TT =Gie)) (4-14) 


has been computed as that portion of the F(S)-insolation incident at 


Bie ssea Surface just prior to absorption by the surface. Note that 





no absorption in air has been included in the computations of Table VII, 
and that the only absorption permitted is that implicit on IS10. The 
final line under each column of Table VII depicts the composite disposi- 
tions computed by means of the weighting scheme of Eq. 4-12 as applicable 
Pormtne case 2/3-CL. 

Table VII. A sample listing of F(S) insolational disposition computed 


at gridpoint (2,2) using equations listed in Sec. IV.B. Insolational 
verdes in ly min -. 


(CL, CL.) Weighting 16S) REFS STRAN 
(0,0) e050 so203 USES, ~4117 

(1,0) - 0000 ~ 1969 meoy 2st 

(Or) - 1609 vlDOL ~ 3342 ree Ae, 

Ci) - 0000 sti0 23742 -1414 
COMPOSITE VALUES 22000 aM Ti MS) 2j BS) 


C. DISPOSITION OF F(A) INSOLATION 
The fractional portion of the solar insolation subject to absorp- 
tion by atmospheric water vapor and carbon dioxide are covered in the 
following subsections. 
ite Clear—-Sky Case (0,0) 
The Manabe-Moller absorptivity function provided the necessary 
absorptivity values for the key layers in this case. The form of this 


iesorptivity function is 


a(2, k) = .271(u(2, k) Sec z]°?> (4-15) 


Here a is the absorptivity applied to the pressure-scaled water vapor 


mass between levels 2 and k (Fig. 2) along the zenith slant-path angle z. 


ae 





The resultant absorbed insolational energy in the partular layer (2,6) 


is then given by the Manabe-Moller relation 


A 


Pie OlN=Oo27 tr (A hu t2,6) Sec zien (4-16) 


The two layers of interest in which absorption was computed were (2,6) 
and (2,10). The absorbed insolation in the layer (6,10) was then com- 


puted by 
A(6,10) = A(2,10) - A(2,6) (4-17) 


Water vapor mass above level 2 was assumed negligible in the F(A) dis- 
position of the solar insolation. 

By subtracting A(2,10) from F(A), the direct transmission of 
F(A) insolation impinging at the earth's surface was determined. The 
transmission of F(A) insolation is then further reduced by the transmis- 
Sivity of (1-a(G)), after surface-reflectance, which leads to the earth- 


absorbed result 


. 303 
} 


TAVO(O, 0) = F(A) 11-.271[0(2,10) Sec 2) (1-0 (G) ) (4-18) 


The transmitted energy impinging upon the earth just prior to absorption 
is 
TRANAO, 0) = 1A10 (0, 0)7 (1=a(G)] (4-19) 
yee Loudy Cases 
In order to compute meaningful dispositions of F(A) insolation 


in cloudy-sky cases, cloud reflectivities and cloud absorptivities 


after C. D. Rodgers (1967) were utilized. The reflectance-values used 
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are different from those suggested by Rodgers for the F(S) wavelengths. 
The cloud reflectivities used here are RA(1l) = .46 and RA(2) = .50. In 
this case there are also cloud-absorptivities to be considered. These 
were taken as A(1l) = .20 and A(2) = .30, respectively. In the following 
discussions the cloud conditions are considered totally overcast in the 
eer Or layers indicated when the notations ((1,1), (0,1), and (1,0)) 
are utilized. 

Schematic representation of the computations performed in the 
wMeetous cases ((1,1), (0,1), (1,0)) are displayed in Fig. 3. This figure 
indicates the theoretical parameters required in computing the disposi- 
tion of incoming F(A)-insolation from level k=2 to the earth's surface 
(k = 10). The equations which relate to the parameters in Fig. 3 and 
the other cloud configurations are listed in Appendix B of Spaeth's 
eiudy (1975). 

The symbols A24, A46, A68, A89, and AQ10, etc., represent the 
insolation absorbed in the layers concerned. Symbols F2, F4, F6, F8 and 
F9, etc., depict the streams of insolation passing through the indicated 
level. A vertical arrow implies the direction of insolation passage, 
i.e., ¥ denotes downward insolation, t+ upward-reflected insolation, and 
¥¥ downward-reflected insolation. The absorption quantity A(6,8)¥¥, for 
example, indicates absorbed energy remaining in (6,8) from a downward 
reflected beam, 

Since multiple reflections occur between the earth's surface 
and a cloud base, or between cloud-layers in this model (Fig. 3), the 
insolational amount reaching the lowermost reflecting surface is 
affected by the downward reflected beams. Additional insolational 


amounts remaining after more than two reflections are negligible, and 
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Figure 3. Schematic representation of F(A)-insolation disposition 


in the case of two overcast layers (cloud-case (1,1)). 
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are not included in the model computations. Note that insolation re- 
flected upward from a lower interface, either cloud-top or ground, to 
the base of a higher cloud deck was not subjected to absorption by or 
transmission through the cloud. This simplification resulted in slight- 
ly reduced F(A) insolation-reflectance to space. 

A final point to note is that there are three contributions to 
the absorption of solar insolation within a layer between any two re- 
flecting surfaces. In Fig. 3, these contributions are identified by 
arrows which indicate the portion of path being crossed. 

For each of the possible overcast-cloud combinations, the 
transmitted F(A) insolation arriving at the earth's surface may be de- 


fined using the following notation (TRANA) as 


TRANA (1,1) = FlLOY + FLOVY (4-20) 
TRANA (1,0) = F1OY + F1OVY (4-21) 
TRANA (0,1) = F10¥ + F1OVY (4-22) 


where the right side parameters of (4-20), (4-21), and (4-22) are 
derived in Appendix B of Spaeth (1975). 
The F(A) insolation absorbed by the earth in each cloud-case 


considered is derived from 


ITA10(1,1) = F1O¥ (1-a(G)) + F10¥Y¥ (4-23) 
IA10(1,0) = F1OY (1-a(G)) + F10VY (4-24) 
TA10(0,1) = F1O¥ (1-a(G)) + F10OVY (4-25) 


In the last three equations, the quantity F1l0+4*¥ is small enough in each 


case, that no further reflections from the earth were considered. 
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Pec Onpestte ein) Layver=Absorptions and Surface- 
Absorption Insolation 


As has been previously discussed, the standard grid-area weight- 
ing scheme of this study was applied to obtain composite values of the 
absorbed F(A) insolation in key layers and also within the earth's 
surface. The weighting factors applied to the corresponding overcast- 
combination absorption quantities provided the following composite 


results: 


m26(CL(1),CL(2) ) A2ZG(ORO)) W(O,0) + A26i0,0) Wl, 0) 


+ 


A26(0,1) W(0,1) 


+ 


AZO re) WAL) (4-26) 


MeO C Cin) 7CiL(2)) =AGLO(O,0) W(0,0) +A610(1,0) W(1,0) 
TROLO (O71) Ww (0,1) 


EAGLOWI, 1) W(t, 1) (4-27) 


meee (1), Cci(2)) =£LAl0(0,9) Ww(0,0) + IALO(1,0) w(1,0) 
+IA10(0,1) W(0,1) 


reo ew lt} (4-28) 


The weighting factors W(0,0),...,W(1,1) were first listed in Eqs. 2-6 
fem. ,d). ALSO A26(0,0), A610(0,0), IA10(0,0) are given in each clear- 
sky case (0,0) about each gridpoint by Eqs. 4-16, 4-17 and 4-18, 
respectively. 

The results for the absorption in layers (2,6), (6,10) and at 
Mae SUrface are shown in Table VIII for the four cloud cases [(0,0), 
wm,O), (0,1), (1,1)) and for the composite case based on the weighting 


maeeors W(O,0),;...,;W(1,1). 
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Table VIII. A sample listing of F(A) disposition parameters computed 
at gridpoint (2,2) using equations listed in Sec. IV.C. Values listed 
ere in ly min Ll. 


(CLs CL.) Weighting A26 A610 AT10 REFA TRANA 
GOO) Ss oa0 ~0509 - 0670 SLL O3 ~03 14 ale is 
le, 0) O20 2052) e032) ~0424 sLO17 SODeS 
(0,1) 2413 OS 2S pS $0211 ~ 0641 20207 
(1,1) OE0) fOaZi 70567 20102 LO? sOl2o 

COMPOSITE VALUES ,0514 POs L - 0888 (0392 el 260 


In the computational scheme indicated by the entries of Table 
VIII, the reflected F(A) insolation to space has been depicted by the 


symbol REFA, and its values follow from 
REFA = F(A) ~- A26 - A610 —- IA1O (4-29) 


whereas the TRANA dispositions are given by Eqs. 4-19,...4-22, respec- 
tively, or by its weighted-mean value in the case of the composites 
Mast line, Table VIII). 

4. Absorptivity (ABA) by Layers 

Here the (fractional) absorptivity as well as the actual insola- 

tion values absorbed in the layers are considered. In the computation 
of absorptivity, which is fractional absorption, the total undepleted 
insolation at the top (k=0) is used. The following equation was utilized 


@aeenis calculation: 


FADJ = 2.00(x/r_)~* COs. Z (4-30) 
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The absorptivity of the troposphere ABA was computed from the ratio of 
the insolation absorbed in the troposphere to the insolation incident 


at the top of the atmosphere rather than at k=2: 


_ A(2,6) + A(6,10) : 
US ee (4—31) 


Do ALBEDO (ALB) OF THE EARTH=TROPOSPHERE SYSTEM 

In considering the planetary albedo, the reflected insolation of 
the earth-troposphere system in both the F(A) and F(S) insolational 
regions must be recalled by the program. Thus REF is computed at each 
gridpoint as the sum of the reflected insolation energy in F(A) (denoted 
REFA in Eq. 4-29) and the reflected part of F(S), previously denoted 


REFS in Eq. 4-13: 


REF = REFS + REFA (4-32) 


Finally the planetary albedo is related to FADJ through 





(4-33) 


E. COMPOSITE ABSORPTIVITY (ABG) BY THE EARTH-SURFACE; COMPOSITE 
ATMOSPHERIC TRANSMISSIVITY (ATRAN) 


ie absorptivity (ABG) of Earth 
By summing the weighted values of F(S) and F(A) portions of the 
incoming insolation entering the earth, the total insolation absorbed 
at the earth's surface was computed. This quantity when divided by the 


extraterrestrial insolation gave the fractional absorptivity (ABG) of the 


De 





earth's surface. The equation for ABG was 


16/21) @ Mis mes EA @) 
ABG = PADI (4-34) 


where IA10, IS10, and FADJ were defined previously by Eqs. 4-28, 4-12, 
and 4-30, respectively. 


Ewe Leansmissivity (ATRAN) of the Troposphere 


Also computed with the total insolational energy TRAN, incident 


Memmeremcalth S Surface just before absorption by the surface. This cal- 


culation is given by 


TRAN = TRANA + STRAN (4=35) 


Here STRAN = (IS10/(1-a(G)) as previously defined in Eq. 4-14, and 
a(G) 1s given in Eq. 4-7. TRANA has also been defined as the weighted 
value of TRANA(0,0), TRANA(1,1), TRANA(1,0), and TRANA(O,1) given by 
ieeweeont9,,,..,22. Note also in justification of STRAN that the four 
Beees tor IS10 of Eqs. 4-8,...,11 each have the common factor (1-a(G)) 
in the numerator and therefore each transmitted F(S)-insolational com- 
ponent available at the earth just before absorption needs only to be 
divided by (1-a(G)). TRAN may thus be viewed as the total insolational 
energy incident at a pyrheliometer located at earth. The (fractional) 


transmissivity of the troposphere (ATRAN) is then computed from 


ATRAN = TRAN/FADJ (4-36) 


Note finally that the major dispositions of the total insolation at 

the indicated map times have now been identified by the fractional 
values, ALB, ABA or ABG, and ATRAN, each of which is a fractional value 
mepresenting the reflectivity (albedo), absorptivity of air or earth, 


Or atmospheric transmissivity as the case may be. 
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5, Computational Check 
The computational scheme utilized in this model was checked con- 
tinuously by summing the fractional values ALB, ABA, and ABG at each 
gridpoint. The value in each case must total .96, Since as previously 
noted the attenuation of solar insolation was taken as four percent as 


it passed through the stratosphere. 


Pee oLATISTICAL ANALYSIS 
In order to substantiate some of the computations performed in this 
section, several of the most important items computed were statistically 
analyzed uSing linear regression computer programs from the BMD set of 
Sracistical programs (Dixon, 1973). 
lee Clear Sky Cases 
Using ALB (0,0,z) which is the clear sky case of albedo, as the 


predictand, and log, Sec z as the predictor, the following best-fit 


10 


equation resulted 


ALB(0,0,Z) = .07856 + .38689 log, )sec Z 


Ry =. VOLG (4-37) 


where Ry represents the multiple correlation coefficient. 
This result was expected over the ocean where a(G) and a(R) 
involved logarithmic dependence on Sec z. The sample-average values 


of ALB and log, ,sec z were 


ALB(0,0) = .1295 Log, psec Z/—~el3L6 
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Other clear-sky regression tests made use of the water-vapor mass path 


length (M) defined by 


{2 


M = (U Sec Ne log, ,(U Sec 2) (4-38) 


This parameterization of water vapor mass is Similar to that developed 
by Hanson (1971), who used a similar M as a predictor in his empirical 
formulations of ABA for both clear and partly cloudy sky cases. In this 
study the two regressions attempted using M as a predictor were with 
ABA(0,0,M) and ATRAN(0,0,M) as predictands. The best-fit equations 


resulting were 


PBA(O,D,M) = [1026 + 2.0281 M 
Ry = 4 Sheva (4-39) 


ATRAN (0,0,M) eo One OS ool 


Il 


R 
M 


Il 


OU J (4-40) 


The absorptivity in the form given by (4-39) agrees closely in form 
with that obtained by Spaeth based on 16 January 1974 data and with 
that of Hanson (1971) based on pyrheliometric data. The means of 


ABA(0,0,M) and ATRAN(0,0,M) for the original 93 gridpoint sample were 


I! 


ABA (0,0,M) 1433 


ATRAN (0,0,M) - 7488 


M =1.44376 (gm pe ae 


2. statistical Relationships Between ALB, ABA and ATRAN 
in the Cloudy and Clear Sky Cases 


In addition to computing the clear-sky value of ALB, it was also 


possible to compute the composite cloud value of ALB. A regression was 
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formed showing the relationship between the ratio ALB(CL(1), CL(2))/ 
fe1O,0) as the predictand and the total opaque cloud cover(CL) as the 
predictor. CL in this model is specified by Eq. 3-12, and gives a good 
approximation of the effective total cloud cover by the two layers of 
cloud amounts CL, and CL., by Eq. 2-4 used in this model. As explained 
eemigeer in this study, a one-third reduction of both CL(l1) and CL(2) 
after initial determination by Eq. 2-4 was also tried in all computa- 
muons involving cloud amounts. This test represented an attempt aie 
tuning the cloud model for radiative calculations, and gave results in 
closer agreement with the latest literature on satellite reflectances. 
Thus all subsequent regressions of ALB(CL)/ALB(0,0) upon CL, computed 
after Eq. 3-12, will include best-fit equations formulated using a CL 
based upon the 2/3-CL model (Eq. 2-5) and the full-CL model (Eq. 2-4). 
A small letter subscript ‘a' in the equation number or symbol 
will indicate those computations uSing 2/3-CL amounts while a small 
letter subscript ‘b' will indicate those utilizing the full-CL model. 
The first regression equations tested were those for RES ME) 
and gue cl) Lor the respective cases of 2/3-CL and full-CL cases. The 
results are given in (4-4la,b). The symbol Ry once again signifies 


the multiple correlation of the statistical regressions. 


ALB(0,0) {1 + 2.9915 CL + .5873 CL“} 


ALB (CL(1)CL (2) ) 


R mes} aia is. (4-41 (a) ) 


M 


ALB (O70) [1 + 227654 cL + .8308 CL“) 


ALB, (CL (1)CL(2)) 


R - 9846 (4-41 (b) ) 


M 
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The sample means of the values in Eqs. 4-4l(a,b) are 


Eq. 4-4la Eq. 4-4lb 

ALB(CL(1) ,CL(2)) = .3163 ALB(CL(1) ,CL(2)) = .3893 
ALB (0,0) = 1295 ALB (0,0). = ,1295 
CL, = 4320 cL = 5890 


The high value of the correlation coefficient in Eq. 4-41 shows strong 
dependence of ALB upon CL in either cloud parameterization case. Also 
note that ALB(CL(1),CL(2)) values computed here are clearly higher than 
reported by Raschke et al (1973), where ALB = .26 for essentially the 
same set of gridpoints during the NIMBUS III, period 16-31 July 1969. 

The 2/3-CL cloud parameterization gives albedo results which are closer 
to Raschke's results, suggesting that the full-CL amounts of Smagorinsky 
(Eq. 2-4) are too high for radiational computations. A second regres- 
Sion was developed between ABA(CL) /ABA(0,0) versus CL. The following 


best-fit equation resulted: 


ABA (0,0) [1+.4667CL - .0148CL"] 


ABA_(CL(1) ,CL(2)) 
BL = .9565 (4-42 (a)) 


2 
ABA (0,0) {1+.5606CL - .1238CL ] 


ABA, (CL(1) ,CL(2)) 


i = ,9597 (4-42 (b) ) 


Mean statistics for these regression cases are 


Eq. 4-42a Eq. 4-42b 

ABA_(CL(1),CL(2)) = .1721 ABA, (CL(1) ,CL(2) = .1827 
ABA (0,0). = 1433 ABA (0,0) = .1433 
cin = .4320 cue = .5890 
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From the means of (4-42), it 1S seen that the model specifies an in- 
crease in atmospheric (solar) absorptivity with increasing cloud cover. 
This result iS in agreement with studies of Warner (1974) and Spaeth 
(1975) and with the earlier study of London (1957). 

An analogous statistical regression was then developed for the 
cloudy-sky transmissivity relative to the clear-sky transmissitivity. 


The resulting regressions were 


ATRAN _(CL(1) ,CL(2)) ATRAN (0,0) [1—.6522CL} 


R 


e .9982 © (4-43 (a) ) 


II 


ATRAN, (CL(1) ,CL(2)) ATRAN (0,0) [1-.6647CL] 


R 


M oe (4-43 (b) ) 


and the corresponding mean values of the transmissivities are listed 


below: 
Eq. 4-43a means Eq. 4-43b means 
ATRAN, (CL(1),CL(2)) = .5389 ATRAN, (CL(1),CL(a)) = .4576 
ATRAN (0,0) = .7488  ATRAN(O,0) = .7488 


The parameter ATRAN(CL"') /ATRAN (0,0) given by (4-43a) is in very good 


agreement with the statistical result of Savino-Angstrom, which was 
QO = Q, (1 ~ .66CL) 


for the tropical ocean areas. The coefficient of CL in the last result 

1S known to be a decreasing function of latitude (Budyko, 1958). 
Similarly, the parameters ALB(CL')/ALB(0,0) and ABA(CL')/ABA (0,0) 

behave consistently with respect to increasing CL, with somewhat more 


realistic mean values for ALB(CL') and ABA(CL') than for the corresponding 
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full-CL parameterization. Hence, in the mean, the statistical analyses 
performed seem to show agreement with observational results of other 
investigators. Comparative computations of the atmospheric transmissi- 
vity, particularly in the 2/3-CL case appear to give better agreement 
with ground-based pyrheliometry. 

A major test still to be made is the comparison of the lati- 
tudinal-distribution of ALB_ (CL) and ALB, (CL) with that of the satellite 


miipeado of Raschke et al (1973). 


fee aubEDO COMPARISONS WITH PUBLISHED RESULTS 

The tropospheric albedo computations of the solar-insolation model 
of this section have been presented for both full-CL and 2/3-CL cases, 
respectively. These albedos have been interpolated to whole multiples 
of 5° latitude between 20S and 65N. These computations were made for 
each of the four oceanic meridians and the resulting albedos averaged 
across the meridians are presented as a function of latitude in Table 
IX for both the full-CL and 2/3-CL cases, where comparison is also made 
with satellite albedos of Raschke et al (1973). For sake of consistency 
the Raschke albedo data have been extracted from atlas presentations 
for the period 16-31 July 1969 to 5° latitude grids, on the same four 
meridians, before applying the zonal-averaging process. The total 
Opaque cloud-cover amounts for the two cases are listed under the 
headings "full-CL" and "2/3-CL" corresponding to the composite cloud- 
Sever formula, Eq. 3-12. 

The primary results to be noted in Table IX are the quasi-constant 
melues Of ALB(RAS) between 20S,...,35N, with individual latitudinal 


values not greatly different from the tropical-subtropical mean of 
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Table IX. Comparison of zonally-averaged albedos as computed by this 
Study for both full-CL and 2/3-CL models for 16 July 1974 with those 
of Raschke et al (1973), and the corresponding values of full-CL and 
@ie2/7S-ClL, each as a function of latitude. The bottom line lists the 
cosine-weighted average of the parameter listed in the tabular column 
above. 


Albedo Albedo Albedo Full-CL 2/3-CL 

Lat. full-CL 2/3-CL Raschke AVG. Avg. 
20S moO 2296 20 SSS ie: «269 
15 2347 7200 ais: ~456 Bac 
10 SAMY ~26e - 16 EAST) 7250 

5 2 2345 ~26 7o16 498 
0 eo 7203 224 -615 ~445 
5 - 404 0324 a Papo 2054 
10 S00 a Pas: ool! - 460 
nls, CESS, SPAe 224 022 - 430 
20 7SO2 oo ay so38 431 
25 7350 2209 eal Seas: oo 
30 ao) eae? ao soul «425 
ope 2925 ABS. 24 ooo -3035 
40 one? Pas, 720 - 546 soed 
45 2205 0341 woe Sleyiall - 566 
50 oO ~309 234 sOo2 a5. 
55 So 463 ae 948 mele) 
60 7526 - 288 -40 2035 2459 
65 413 aoe3 noo +862 A SNS 0) 
Wt.Avg. .364 oO Ceo soe ~405 
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ALB (RAS) = .22. For most of this latitude range the computed values of 
ALB (2/3-CL) exceed those observed by Raschke (Table IX) by the frac- 
tional amounts .05 to .10. From 40N poleward ALB(RAS) and ALB(2/3-CL) 
are in fairly close agreement, with individual values close to the local 
average of 0.35. Note that at latitudes 60, 65N, the ALB(RAS) results 
are slightly greater than ALB(2/3-CL) which may be due to an effect 
associated with the large mean zenith angle at these latitudes. 

As noted earlier in Table V, at latitude 55N there is a localized 
maximum of CL‘ associated with the largest computed value of ALB(2/3-CL) 
= .463. The later value is not verified by the two-week average obser- 
vation of ALB(RAS) = .37. Overall, the computed values of ALB(CL) 
exceed those of ALB(RAS), latitude-by-latitude, in much the same degree 
as indicated by cosine-weighted means of Table IX, which are also 


listed below: 


ALB (full-CL) = .36 
ALB(2/3-CL) = .30 
ALB (RAS) = .25 


To summarize, there 1s a well-defined preference for the 2/3-CL 
over the full-CL model in giving albedos closer to those observed by 
satellite for the same general period (16-31 July 1969) as reported by 
Raschke et al (1973). An area of difficulty even with the 2/3-CL 
model occurs in tropical and subtropical latitudes where Von der Haar 
and Hanson (1969) have reported the reality of smaller than suspected 
cloud-covers and of smaller values of the resultant global-albedo. At 


these latitudes (tropical and subtropical) it seems doubtful that the 
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large-scale CL'-parameterization, following Smagorinsky, is as applicable 
as in higher latitudes. Rather some method of specifying the distribu- 
tion (at gridpoints) of isolated cumulus towers may be appropriate. 
Finally, it 1s questionable whether vertically-structured convective 
cloud elements will have the high reflective capability attributed here 
to clouds which are depicted as existing in horizontal layers with a con- 
peemeererlectivity 2.5 < R< .65 as specified in the radiation model. It 
may well be that the vertical shafts between convective elements may be 

a more efficient focus of solar energy to the earth than is perceived by 
ascribing a constant cloud albedo to the cloud-layer. An alternative 
approach would be to reduce the two-cloud layering in the tropics to a 
Single layer of clouds. 

At any rate with the two-layered cloud-model considered in this 
study, the only feasible selection from the results deduced here is the 
2/3-CL parameterization. This selection is carried over into the compu- 
tations of Section VI. A further test of the cloud-parameterization 
effect on the computation of net radiation of the earth-atmosphere 
system, by comparison with that observed by satellite, is made in Sec- 


mom VIL. D. 
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V. SENSIBLE AND LATENT HEAT TRANSPORT 
AQ STHE orn th INTE REACE 

fi. GENERAL PURPOSE 

The model used here to describe the turbulent transports of sensible 
and latent heat across the sea-air interface was basically that already 
in use in the operational FNWC primitive equation model. This model has 
been discussed in detail by Kesel and Winninghoff (1972) and by Kaitala 
(1974). The primary purpose of the adaptation of the turbulent flux 
model in this study was to test the magnitudes and directions of these 
transfer rates in comparison with those of the radiational transfer model 
which has been discussed here in Sections III and IV. The final purpose 
of the inclusion of the turbulent flux model was to conduct a diagnostic 
computation of the heat balance within the oceanic and atmospheric 
layers for 16 July 1974 under the combined operation of the radiational 


and turbulent flux models. 


B. WiINDSPEED COMPUTATION IN THE TURBULENT FLUX MODEL 
At all of the 93 gridpoints tested, the geostrophic windspeed ie was 


computed by use of the standard formula 


_g _ Az 
we Zs Auy ae 
; AT 
with —— (36, 400° Sin 9 (5-2) 
a 1.86603 
poli oO 
/ 2 A 
= _ 4- - = 
and NZ (2, 2.) (2. Z,) (5-3) 
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In Eq. 5-1, da = 381.0 km is the nominal FNWC grid-spacing (true at 60N) 
and (ES) is the true-earth spacing used in the centered-difference com- 
Pueati1ons Of the contour gradient (cf., Fig. 4). 

In (5-1) and (5-3) AZ is the contour-difference, centered in most 
cases on the gridpoint under consideration. However, at the map edge 
gridpoints of the meridians, it was only possible to compute a forward- 
differenced version of Le ora Pig. Ae, ens at edge=point "0", Eq. 
(5-4) below = used. The spacing . in Eq. 5-4 is then taken as one- 


half that in (5-1), so that 





_g AZ 7 
Vg Se ayn) (5-4) 
with 
— = 2 e> 2 , = 
AZ = V(Z, zo) se (Z. Zo) (5-5) 


All geostrophic windspeeds were then converted to surface windspeeds us- 


ing the following empiricism due to Langlois and Kwok (1969) 
V= .8V + 2.2 (5-6) 
Ss g 


where ve = surface windspeed (mps) and We is the 1000 mb geostrophic 
windspeed. The factor 2.2 in (5-6) is an empiricism introduced to in- 
clude the effects of gustiness upon the measured (steady-state) wind- 
speed. 

Pimetemintervea et ho degrees of latitude, the value of the 
Bertolic parameter was arbitrarily set at f = + .25 x MOmccce follow- 
mag Kaitala (1974), while for | > | Palo je eene actual value of £ for the 


particular latitude was used. 


AN 





enw 





Figure 4. Section of FNWC polar stereographic meridian illustrating 
the method of obtaining contour gradients in the vicinity of a grid- 
point. The values of contour-heights correspond to the 1000 mb D- 
values in the neighborhood of the gridpoint. 
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SOE NSIBLE HEAT FRANSPORT 
The details of the adaptation of the FNWC model for sensible heat 


transport H, utilized in this study follow the description of Spaeth 


ia 


(1975). The formulation of H, iS given below in terms of bulk para- 


If 
meters of the surface and of the planetary boundary layer (10,8) of 


Pg, 2: 
Hoe Owe Ce eae - TT) 
-j 
ee Si ey See a 
-1 -1 
where oe =,239 cal om (deg. K) 
3 
and oe = 1.4 x 10 (after Weiler and Burling, 1967) 


The equation utilized in this study to compute i may be shown to be 


(Spaeth, 1975) 


2K O5 7” 
Covi eT + -2kKy 
T z Ss og - AZ @. (5-8) 
Zs 
—— + 
NZ, Cp we 
where 
MOCO. 236 
oe = To S00? (5-9) 
* 
by K 
NS (5-10) 
on G 
1 + a* (> 3 J) 
oe ae 
R 
qd - 10 
KZ = Bo 250 = Ge een 9? (5-11) 
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* * 
In 5-7(c), the turbulent parameters K anda have been determined 


empirically by FNWC 


k = 
K = qowenaeee 
* = 
Bl 85) x Oem) 1 y 
ue - T in Eq. 5-7 is the temperature lapse for the thin layer between 
x 


eieesurface and the top of the constant flux layer (cf., Fig. 5). The 
formula for Ty deduced here required a constant value of sensible heat 
Hp within the constant flux layer (CFL) at each gridpoint with Hp then 
decreasing linearly with pressure to Hp = 0 at k=8, and then remaining 
zero above k=8. This linear decrease of Hy with pressure is compatible 
with the condition of constant convergence of the sensible heat flux per 
gram in the turbulent boundary layer (x,8). Thus, the constant value 

of the convergence of sensible heat as given by Eq. 5-7 within the layer 
(8,10) may then be considered to be applicable at level k=9 as depicted 
marie. O(a), after Spaeth (1975). 

After computing ui Dy che D—G, Hp was then computed at each grid- 

point using the bulk transfer equation 5-7. The gridpoint values of H 


r 


are used later in calculations of heat balance in Sections VI and VIII. 


D. EVAPORATIVE HEAT TRANSPORT 

Similar to the case of sensible heat convergence, the latent heat 
transport by turbulence was assumed to be subject to constant evapora- 
tive flux convergence throughout the layer 800-1000 mb. The total 
amount of latent heat removed by evaporation from the ocean surface was 
essentially as modeled by Kaitala (1974) using bulk transfer theory. 


It should be recognized, however, that the condensation and consequent 
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mealized latent heat may occur at arbitrary levels, i1.e., in general at 
levels higher than k=8. Thus, in following the FNWC model, we have 
not introduced E as a heating rate at level k=9. 

The basis of the turbulent latent heat model adopted is shown in 


Fig. 5(b), where 


E (900) = 1/2 E (Bulk) (o-) 


and where the surface layer bulk transfer of latent heat is given by 


E = E(Bulk) = Le Cv. (q_-a,) (5-13) 


Evaporation at level k=9 may also be computed using 


KGa Cj) 

= Xx 9 

E = i K (5-14) 
Ae ye easy: 


ww 


where K is given in Eq. 5-10. In (5-13) and (5-14), Lis the latent 


heat of vaporization. 


Pea l6.75 )— 7.60. Tho 


alg! Pelicm) ies.) seo 14 and 5=—15 do and qd, are the mixing 
ratios at the levels k=9 and at x, respectively, whereas qe is the 
saturation mixing ratio at the surface temperature ae The tactor 1/2 
which appears in the right side of (5-12) reemphasizes the concept of 
Sootant latent-heat convergence in the layer (x,8), (cf£., Fig. 5b)). 
The solution for qd, 1s obtained by combining Eqs. (5-12), (5-13), 


ame (5-14) in the form 


~ 
? fe 


K 
Ss ve Ti0 Pre Az to 
Ll a a (S215) 
+ ——— 
CD we ¢ he 


- 7 . 


Wienre  /\Z = Zo 10 
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Rep licderom Of BGs.) 5-13 and 5=15 makes possible the computation of 
E at each gridpoint. As in the case for T, the factor of Pee) el 0 
5-15 was taken to conform to the value K as uSed by Kaitala (1974) so 
that the turbulent transport results obtained here will be compatible 
with those computed operationally by the FNWC primitive equation model. 
As already noted, this was accomplished in this study by using a value 
Seelk* = 2 xX oemace cme aS equivalent to the K*¥-value of Kaitala (1974). 
Thus computation of E by Eqs. 5-13 and 5-15 was equivalent to that using 
Kaitala's scheme (1974), the only change being that the vertical distri- 
pueden Of & Implicit in Fig. 5(b) is slightly different from that of 


the constancy of E everywhere in the planetary boundary layer, as postu- 


Matted by Kaitala (1974). 


E. LARGE-SCALE TURBULENT HEAT FLUX ACROSS THE SEA-AIR INTERFACE 


The values of H, represent the sensible heating rates for the atmos- 


Ni 


pheric layer (8,10) and were computed for each gridpoint along the four 


meridians. The values of H, in general turned out to be large and posi- 


IP 


tive in the Southern Hemisphere reflecting the computed positive values 


was predominantly negative 


r 


(ere, Figs. /,...,10) and therefore contributes negatively to the net 


of ue - Ty - In the Northern Hemisphere, H 


atmospheric heat budget in the layer (8,10). 


The latent heat term (E) was large and positive relative to H, in 


if 
the Southern Hemisphere and in the latitude band 0-25N. Regardless of 


rT combined values of E + Hp seemed 


to give reasonable heating rates in the atmosphere in these latitudes. 


localized variations in the sign of H 


On the other hand, the combined E + H, acts as heat-loss contribution 


i 


to the troposphere in Northern Hemisphere latitudes ¢ > 30N and is of 


ia 





the same sign as the sensible heat~loss contribution (in these lati-~ 
tudes) applied to the layer (8,10). 
In those latitudes (Southern Hemisphere and the Trade Wind Belt) 


where E + H,, acts aS a positive contribution to the heat budget of the 


r 


tropospheric column, it 1S not specified in this study precisely where 
the latent heat source E will be realized. The use of the equations of 
conservation of mass, momentum and of water-vapor mass, would resolve 
the latter question, but this was considered beyond the scope of the 
present work. 

mie turbulent flux computations of Section VI for Hy and E just out- 
lined represent minor modifications of the FNWC turbulent boundary 
layer model which has been appended to the major concept of this study 
namely the radiative model. The purpose for inclusion of the E, Hy 
computations has been primarily for consistency tests, that is, to 
determine whether net radiative flux at the surface is of the proper 


magnitude, relative to E + Hp « 
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Vig MEREEPTONAL CROSS-SECTIONAL DEPICTION 
OF THE HEAT-BALANCE COMPUTATIONS 

A. GENERAL 

The general design of this section is to utilize all of the computa- 
tional concepts discussed in Sections III, IV, and V in the computations 
for a single time-step in the heating model developed for use in the FNWC 
prediction model. After testing FF2 and ALB with the corresponding 
values of Raschke et al (1973), it was decided that only the computations 
by the 2/3-CL parameterization would be displayed in the meridional 
cross-sections (Figs. 6,7,8,9,10). The appropriate calculations were 
performed at each gridpoint of the four meridians used for presentation 
purposes in this section. However, the same type of computations using 
the full-CL parameterization were also made, but not presented in this 


Section. 


B. GEOGRAPHICAL REPRESENTATION OF THE HEAT-BALANCE DISTRIBUTION 

The FNWC gridpoint processed analyses for OQOOOGMT, 16 July 1974 were 
used at the three Pacific cross-sections, while that for 1200GMT, 16 
July 1974, were used for the single Atlantic meridian. This was done 
so that the set of gridpoints was considered to be subject to the actual 
radiative-transfer calculations involved for these specific times. 
Figure 6 depicts in symbolic language the key to the computational en- 
tries in Figs. 7,8,9 and 10. This symbolic list presents initially the 
computations made at each gridpoint (I,J) radiative sounding data in 
the form of Table I(b). The computations proceed from the top of the 


troposphere to the ocean surface and include the sensible and latent 
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heat transfers, by the means of the computational procedure of Section V. 
For purposes of climatological data comparison, Figs. 7,8,9 and 10 were 
developed by interpolating gridpoint results to whole 5° latitude incre- 
ments. The interpolation routine to this gridpoint spacing made use of 


the Lagrangian cubic interpolation scheme (after Spaeth, 1975). 


(f= 2) Glas eebe 4) (I-1) (1-3) (1-4) 


QO(I) = Q, MED les) eR Q., (2-1) (2-3) (2-4) 


(oes) 
Cla) aa al) ia ie) hl = 3) 


23 (31) 32) (e=4) 24 (4-1) (4-2) (4-3) 
Finally for ease in reconciling the magnitudes of all heat-transfer 
rates, the time-dependent solar disposition rates have been averaged to 


24-hourly rates. 


C. EXPLANATION OF SYMBOLIC TERMS 
1. Cross-Section at Level k=2 (Fig. 6) 

The discussion of all insolation parameters discussed previously 
in Section IV dealt with the specific time of day that corresponded to 
the hour angle h for the instantaneous time t under consideration. The 
incident solar insolation dealt with is then 

F(S) =S (=)? Cos ez (652) 
m 
In order to avoid reference to specific map times t, the instantaneous 
Solar hour-angles were h = 55°, 10°, 35°, and 35°, respectively for 
Sross—sections 1,2,3 and 4, as depicted in Figs. 7,8,9 and 10. 
QAVE represents the 24-hour average of F(2) and appears as the 


first input symbol in Fig. 6. Its value is considered to be more 
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representative climatologically for the data day under consideration 
iemam & (2) < 


QOAVE is derived by the formula 








Cos Zz 











ORVE = F(2) = — (6-3) 
where 
Cos z = [HSind Sind + Cosd Cosdé SinH]/7T (6-4) 
and 
= arcCos[- tang tand] (6-5) 


Hesewo = 265° is the solar declination for 16 July 1974, and H is 





the hour angle at local sunset at latitude 0. Cos z in (6-3) 1s equal 
to the 24-hour average cosine of the zenith angle, Eq. 4-2. The 24- 
hour time averaging sysctl’ yielding QAVE gives heating results consis- 
tent 1n magnitude with the terrestrial-flux divergences, which change 
miveslighntly with time of day. This conversion to expected daily- 
averaged solar disposition guantities 1S compatible with the determina- 
tion of a heat budget for the een date (16 July 1974). 


Other parameters needed for level k=2 are 








CCseZz 
OREF = REF (t) get (6-6) 


where 


REP Cj —sbe2)) -A126 -A6lO0 —(IA1LOFIS10) (6-7) 


REF(t) 1S the instantaneous solar reflected insolation at a gridpoint 


and QREF is its 24-hour average, assuming that the planetary albedo 
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remains constant for the 24-hour period, although it was calculated for 
the indicated solar time t. This assumption requires that the cloud 
amounts computed at the indicated synoptic times are representative of 
the entire day. 

The same principle will be used with regard to all other solar 
parameters in the conversion from time-dependent values at solar time t 
to 24-hour averaged values. Superior bars ( ) are not used in present- 
ing the averaged values shown in the cross-sectional key, Fig. 6. The 


24-hour average earth-troposphere system balance, BALT, is computed from 


BALT = QAVE - (OQREF + PF.) (6-8) 


for the tropopause level k=2 at the indicated latitude. Net terrestrial 


* 
fluxes, such as F, , were considered to be constant throughout the 24- 


2 
hour period, a valid assumption if the cloud cover remains quasi-con- 
Stante for the period. 
2. Cross-Section in Layer (2,6) 
The following definitions apply for the layer (2,6) as referred 
to in Fig. 6. All the heat transfers shown in this layer are assumed 


to be of radiative character only, as in the corresponding version of 


the FNWC heating model. The averaged radiative cooling rate is given by 


BAL 26 = 026 - F26 (6-9) 
where F26 = terrestrial cooling rate in 
ehitene 1 AS) 
and 026 = daily solar absorption in layer (2,6) 


ws defined relative to A26(t) by a cosine transformation similar to 


EG. 6-3. 
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Beep etoos-ocetlon in bayer (6,10) 


The 24-hour average radiative cooling in the layer (6,10) is 


given by 
BAL 610 = Q610 - F610 (6-10) 


where Q610 is the sum of 068 and Q810. BAL 68 has been taken as one- 
half of BAL610. However, the heat balance in the layer (8,10) has been 


hiodafled for inclusion of Hpe as given by Eq. 5-7. 


BAL810 = 0810 - F810 + Hp = BAL68 + Hp (6— 1s) 


Here Hp is the sensible heat-flux convergence in the layer (8,10). 
F610 is the sum of F68 and F810, the terrestrial flux losses for the 
respective layers. 

4. Cross-Section at Air-Sea Interface (k=10) 


The heat balance at the earth's surface (BALB) is as defined 


in the following equation 


x 
BALB = QABG - F - (H 


We + EB) (6-12) 


r 


QABG is the 24-hour average insolation absorbed by the surface defined 


by 





OABG = QABG(t) (Cos z/Cos z) (6-13) 


Hp and E were computed, respectively, from Eqs. 5-7 and 5-13 and con- 


Sidered constant for the computational day, 16 July 1974. 
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Pee RE DLONAL CROSS=SECTIONS OF THE VERTICAL HEAT BALANCE 
FOR 16 JULY 1974 


Figs. 7,8,9 and 10 as previously explained represent the single time 
step of heating computations for each of the four meridians used in this 
Study. The four figures have been divided into (a) tropical results and 
(b) mid-to-high-latitude results. While the results depicted in these 
cross-sections are exhibited as representing daily-averaged values, they 
are actually based upon heat-computations at the specific map times of 
OOOOGMT and 1200GMT on 16 July 1974. Therefore, for these results to 
be meaningful aS a Stepwise part of the heat package subroutine of FNWC, 
the solar radiative absorption and reflectance terms would have to be 


recoverable as a function of GMT, l.e., 





J QAVE (Cos z/Cos z) (6-14) 





REF (t) OREF (Cos z/Cos 2) (6-15) 


etc. Thus solar disposition terms may then be utilized in connection 
with the one-hour time stepwise application of the thermodynamic equa- 
tion of the set of primitive equations used in the FNWC prediction 
process, for the 2/3-CL parameterization used throughout this section. 
The full-CL heat~-balance cross-sections were also computed, but have 
not been presented here for the sake of brevity. However, certain key 
radiative heating (cooling) rates have been summarized in zonally- 
averaged form for both cloud parameterizations for the purposes of pre- 
senting the mean latitudinal comparisons (Section VII) of the purely 
radiative heat budget of the FNWC data as it was applicable to 16 July 


1974. 
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REECULtO 21g. 6 Orne. 


125°W Longitudinal cross-section for tropical latitudes. 


Falguve 74a). 


Values computed from data for 16 July 1974 for 2/3-CL case only. 
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Refer to Fig. 6 for key. 


125°W Longitudinal cross-section for higher latitudes. 


Figure 7(b). 
Values computed from data for 16 July 1974 for 2/3-CL case only. 
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170°W Longitudinal cross-section for tropical latitudes. 


Values computed from data for 16 July 1974 for 2/3-CL case only. 


Refer to Fig. 6 for key. 


Figure 8(a). 
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Refer to Fig. 6 for key. 


170°W Longitudinal cross-section for higher latitudes. 


Figures (b>). 


Values computed from data for 16 July 1974 for 2/3-CL case only. 
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Refer to Fig. 6 for key. 


145°E Longitudinal cross-section for tropical latitudes. 


Figure 9(a). 


Values computed from data for 16 July 1974 for 2/3-CL case only. 
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Refer to Fig. 6 for key. 


145°E Longitudinal cross-section for higher latitudes. 


Figure 9(b). 


Values computed from data for 16 July 1974 for 2/3-CL case only. 
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Figure 10(a). 
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35°W Longitudinal cross-section for tropical latitudes. 


Values computed from data for 16 July 1974 for 2/3-CL case only. 
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Refer to Fig. 6 for key. 


35°W Longitudinal cross-section for higher latitudes. 


Values computed from data for 16 July 1974 for 2/3-CL case only. 





Vitti ba ODiNAl DISTRIBUTION OF RADIATIONAL 
BPA wer era wor ih, OCREAN-ATMOSPHERE SYSTEM 
A. GENERAL 

The latitudinally distributed results of the purely radiative contri- 
butions, over the ocean-atmosphere system, are presented in Figs. 12(a,b) 
for the 2/3-CL and 13(a,b) for the full-CL cases, respectively. These 
cross-sections show the results after averaging over the four meridians 
considered in this study. The results are displayed in the format of 
Fig. ll. Radiative transfers 068 and 0810 have been combined with F68 
and F810 of Fig. 6 to compute layer-average radiational warming rates in 
the layer (6,10). The turbulent transfer heating rates E + Hp to the 
atmosphere (and the corresponding cooling rates in the ocean) which were 
included in the results of Figs. 7,8,9, and 10 are now omitted in the 
radiative composites presented in Figs. 12 and 13. 

In obtaining latitudinally-distributed means of radiative heating 
rates at each five-degree multiple of latitude ¢$, all values of each 
radiative parameter listed in Figs. 12 and 13 at latitude ¢ in the range 
20S,..-,05N (5 degree increments) were simply averaged over the four 
meridians providing the zonally-averaged value at ¢. At >? = 65N there 
was only one contribution to the latitudinal average, while in the South- 
ern Hemisphere latitudes 20, 15 and 10S, only two values (on XX = 125W 
and A = 35W) of each parameter contributed to the means listed in Figs. 


12 and 13. At 5S, there were three sets of radiative parameters enter- 





ing into 0(¢). Otherwise, there were four values of each radiative 
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parameter entering into the computed radiative mean values Q(>) of the 
cross-sections. Therefore, near the northern and southern boundaries 
of Figs. 12 and 13, the listed values may not be equally representative, 
although the general trend should be reliable. Fig. 14, uSing data 
extracted from Figs. 12 and 13, tends to confirm this statement. 

Figs. 12, 13 and 14 are also useful in determining the relative 


merits of the 2/3-CL model as compared to the full-CL model. 


B. EARTH-TROPOSPHERE SYSTEM RADIATIONAL BALANCE SUMMARY 

The latitudinally averaged distributions of the three key radiative 
Parameters, Ros R, and R as functions of latitude for both the 2/3-CL 
and full-CL cases are shown in Fig. 14. The radiative-balance parameters 
Bowe defined {after Malkus, 1962) as R.(9), the mean radiative-energy 
transfer rate across the top of the ocean-tropospheric system at k=2 
ieererred to as BALT in Fig. 11); Ri(?), the meridional mean radiative 
cooling rate in the troposphere (BAL26 + BAL610 in Fig. 11); and R (9), 
the mean radiative warming (cooling) rate at the earth surface (BAL10O 


in Fig. 11). The relationship between these radiative parameters is 


R = R_ +R (7-1) 


The 16 July portion of Tables X, XI and XII list the zonally-averaged 


values of Ros R and R, as a function of latitude $9 for the 2/3-CL and 
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the full-CL cases. Weighted averages for any Q-value are then computed 


using the following cosine weighting scheme (Eq. 7-2, Spaeth, 1975) 





18 4 
k oe 
i ji 

p2 ; : 2. a ine Pi 

ex — a qm 
wt. AVG). 18 2 wae) 

= Cos 9¢ 

i=1 


where aA = value on meridian j at latitude 5 and kK. is the number k=l, 
-e-e,4 Of meridional observations available for the meridional averages 
ite >.- Note that i=l,...,18 corresponds to the 18 latitudes o.= -20., 
eoe,09. degrees. Tables X, XI, and XII also depict the comparative cross- 
seasonal values of Ros R, and Ro by latitudes for both 2/3-CL and full-CL 
Parameterizations for 16 April 1974, 16 July 1974, and 16 October 1973, 
respectively. 

In Fig. 14, the net-flux terms across the interfaces k=2 and k=10 
are R_ (0) and R(¢?). At k=2 net-flux values obtained from the 2/3-CL 
model (denoted Ro) exceed those by the full-CL model (R_.) at ali tatty 


tudes. There is a localized maximum value of both R. and ie ino tae 


1 2 


latitude zone 40 to 50N, with localized minima at latitudes 45 and 55N. 
Table IX indicates the existence of maximum CL at these two latitudes 
(45 and 55N) suggestive of polar frontal-zone cloudiness and higher than 
normal global albedo. The peak hemispheric values of Roa and Ro» then 
@eeur finally at latitudes 65N in this study. If computational data 
extended further poleward, the overall maximum Ro could possibly occur 
at latitudes higher than 65N, in view of the high declination of the 

sun on 16 July 1974. The same general pattern in the zonal mean radia- 


tive balance R(¢) at sea-level exists for both the 2/3-CL and full-CL 
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All values listed are daily averages in ly min71! and are 
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Zonally-averaged radiational cross-section for mid to high latitudes with 


Essen. 2 (5): 


All values listed are daily averages in ly min72 


RetexytOurig.- ll ter key. 


2/3-CL model. 


and are from data for 16 July 1974. 
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Similar to Fig. ll(a) except full-CL model (after Smagorinsky, Eq. 2-4) is used 


11 radiative calculations. 
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Similar to Fig. 11(b) except full-CL model (after Smagorinsky, Eq. 2-4) is used 


in all radiative calculations. 


Figure ss). 
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Figure 14, Radiational balance at the tropopause (Ro), at the 
ocean surface (R), and in the tropospheric column (R_). Solid 
lines denote computations made with the 2/3-CL (subscripted "1") 
cloud model and dashed lines correspond to the full-CL model, 


(subseriprecam2 ). All computations for for léeduly 1974 in ly 
min-l x 102, 
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cases. However, the sea-level maximum in the 40-50N latitude zone is 
more pronounced than that at latitude 65N. It is to be noted that R, 


for the 2/3-CL case is greater than R, over the entire range of lati- 


2 

tudes under consideration in Fig. 14. 
In particular Fig. 14 for 16 July 1974, shows that the quasi-con- 

stant difference Roy - Ro»5 from the bottom line of Table X(b) is 

-0346 ly “io Likewise, Fig. 14 shows that the net-flux difference 

R, (9) - R, (9) is closely approximated by .0303 ly ain, the result of 

the cosine-weighted difference from Table XI(b). 


From the flux-divergence equation 7-] pertaining to the radiative 


cooling of the tropospheric columns 


the values of Ri and R.2 turn out to be nearly equal, regardless of 
cloud configuration. From the considerations of the preceding para- 


graph, the coSine-weighted mean difference Rai - Ri» is 


—S$> i—-  ——- SO ——,— Or —-l— 2 
et ee ee BS) 0043 ty min 


which also results from Table XII(b). This is a relatively small differ- 
ence conSidering the mean cooling rates involved. Fig. 14 makes this 


same result evident upon examination of the curves of Ri and R2° 


There are consistent features in both the Roi and Ri curves, such as 


the minimum values in the latitude zones -10, to the equator and 20N to 


30N, which make the general results on R, and Ro appear credible. 


I 2 
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Ce CROSS — SEASONAL EFFECTS 

Tables X and XI show the three cross-Seasonal dates affecting R. 
and R centered on 16 July 1974. It is clear that Rai is a relative 
cross-seasonal minimum at all ¢ < 5N for 16 July 1974. A similar 


effect applies to R. thrugh the same time sequence. Intra-seasonal 


2 
maximum results apply poleward of 1ON to the cross-seasonal compari- 
sons of R at the earth's surface in either cloud-case,. 

Boia eomdegmees sOllatitnde in Naples X and Xl, there is a 


systematic trend toward minimal values of both R. and R, at 16 July 


1 


1974. The same effect applies to Roo and R, Onstils date, For this 
limited geographical range near the Southern Hemisphere border, the 
16 October ae having a solar declination 6 = -8.5° latitude corres- 
ponds to the largest positive values of both R. and Rk Tiny elthner ry cloud 
configuration. These cross-seasonal results appear to afford quali- 
tative support to the radiational model tested here. 
Poe COMPARISONS OF NET FLUX AT THE TOP OF THE MODEL ATMOSPHERE 

With SATELULTE OBSERVATIONS 

Based upon observations of Raschke et al (1973) of F2(RAS) and of 
the albedo ALB(RAS) over the latitude-range of this study, the net 


flux parameter at the top of the atmosphere was defined for the NIMBUS 


III data-period (16-31 July 1969) at k=0 


Ee (1-ALB (RAS)) - F2 (RAS) (T-3} 


= x 
RN (RAS) = QAVE* +55 


Raschke's computed values of the RN for this data period used a slightly 
different solar constant than that employed in this study (S = 2.00 ly 


min7ty , Hence for purposes of consistency in comparisons, our average 
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Extbacerrcstrral insolation, OAVE*(2.00/1.92) was used in the defini- 
ron Of SRN(RAS) in Eq. /=3. 

Values of RN(RAS) were computed after Raschke at each five-degree 
latitude intersection on the four meridians where comparisons with the 
present model computations were possible. For the corresponding radia- 
tive model computations of this study, the following analog to Eq. 7-3 


was adopted. 


RNMOD = QAVE* a (1-ALBMOD) - FF2 (7-4) 

FF2 represents a parameter analogous to F2(RAS), values of which 
were compared in Section III.E. on a latitudinal basis. Likewise, 
albedo values were compared on a latitudinal basis in Section IV.G. 

The computational values of RNMOD for 16 July 1974 were computed 
uSing both the 2/3-CL and the full-CL cloud parameterizations. After 
obtaining averages across the four meridians, the latitudinal compari- 
sons are depicted in Table XIII. 

Comparison of RN(RAS) with RNMOD in Table XIII shows that values 


of RNMOD are somewhat smaller than the values of RN(RAS) at all lati- 


tudes > < 45N. The positive difference 


RN(RAS) - RNMOD > < 45N 


occurs because of the consistently small albedos of the Raschke 


observations. This difference is considerably larger for the full-CL 


case, when the model albedos were excessively high. 
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R -values for the earth-troposphere have been extracted from Table 
S 


X and appear in columns four and five of Table XIII. Recall that 


* 


R_ = QAVE (1-ALBMOD) - F, (7-5) 


maecomparing Eqs. 7-4 and 7-5, it 1s clear that the incoming insolation 
used in the computation of R_ has been reduced by 4% of an amount which 


is approximately 





a =i] 
QAVE - QREF = .4148 ly min (7=G) 


when cosine-averaged over all latitudes considered in this study. 


Table XIII indicates that 


— yates 
Ro RNNMOD= 2/73) = 9.025) 8y Min 


holds for the 2/3-CL case. This difference may be reconciled with the 


net-flux differences of Eqs. 7-5 and 7-4, through the differenced 


equation 
R., ~ RN(MOD-2/3) = (FF2 - F,*) - .04(QAVE - QREF) LD 
which reguires that 
(FF2 - ae = .025 + .04(QAVE - QREF) (3) 


; eel 
in ly min for the 2/3-CL case. 


The result of Eq. 7-8 shows the average IR-flux divergence across 


* 


the layer k=2 to k=0 by the F, 


model of Eq. 3-5. The water-vapor and 
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co. absorber masses of the stratosphere are therefore not quite in radia- 
tive equilibrium, because these same absorber masses do not have a com- 
pensating solar absorption effect in the present model. The 4% solar 
absorption by ozone and oxygen are not sufficient to compensate for the 


added water vapor-CO. flux-divergence in the terrestrial waves. 


Z 
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(a) 


16 April 1974 


Mid-Seasonal Dates 


16 July 1974 


(e} 


16 October 1973 


ae ehh Ro2 ol Ye2 ei Neo 
20S -,0229 -.0469 -,1012 ee llza ~0764 -0384 
i> =, OL -.0418 = Otel i -.1010 207385 -0400 
10 -.0002 -.0358 -.0546 =O 9 2 -OsSUs -0430 
5 0171 -.0235 -.0468 -.0656 - 0814 - 0480 
0 -O03050F 2.0079 Aoal sss! -.0154 A (OhSs S78} 5 0)8.31// 
5 BOzGs 20243 sO207 = 00s. sOSLe SOsce 
10 Sol behs: -.0364 20D2E ~0145 - 0484 -0164 
12 OES, =. 0355 ,0052 S02 I 20170 -.0030 
20 -0204 -.0368 - 0627 2UL36 —,0012 =,0322 
25 =,0256—5 ~.0126 - 0825 - 0380 -.0281 -.0588 
30 20047 sU293 ~0945 fOso0 =.09/3 -.0840 
35 - 0488 Os OU ol OS -0607 -.0743 -.0916 
40 -0540 0270 ©1254 .0888 -.0916 -.1036 
45 ~0427 - 0096 «0892 soleul 2! =o od =. 1212 
50 = 010 2S =, 03072 ~1461 oA a)eys: -.1347 -.1395 
>> -0016 =, 0146 eLO63 - 0766 =i SOK =. 1302 
60 -.0428 = Weil coos eplke al2: el oe -,1493 
65 -.0089 —.03/4 Pele)? 0 - 1640 -.1169 -.1096 
Wt.Avg. e002? -,0153 7Oa06 SOLES) 0, —sOUGL -.0320 
Table X. Cross-seasonal comparison of zonally-averaged values of R 


as a function of latitude. The averaged radiational values for the” 
P7o=ci, cloud parameterization are denoted with a subscripted 1 and 
fie f£u11-CL values by a subscripted 2. (Ali values in ly min@}). 
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Mid-Seasonal Dates 


(a) (b) (c) 


L6rApyer) 1974 16 July 1974 Gee rober 1974 

LAT. Tl R eS R, Ry SS 
20S 20s13 ~0514 .05 95 ~ 0547 ~1742 51329 
5 aplonls: SOlsreul 20954 ~0852 ~ 1845 - 1438 
10 e299 Ooo el 7 ~1134 aplsMla! 01524 
5 435 - 1049 shoo? Pl Bessy) BL ZoG faa 2 
0 ~1164 eG ~1741 2508 5 ny, Sle? 
Ss aicoe ~1294 Pal bg a - 1464 S208 7 skuZ2 
10 elo ~ 1460 ~ 2042 clic. 21386 sooo 
5 Aaneyats ~ 1470 S226 - 1866 olG57 ~1465 
20 aZ018 IESE, o2ou3 21959 Pasalal wel OS 
25 22079 LIES 2500 e2225 - 1490 28282 
30 s2053 ee al 7 ae OG S226 ehh 9 70958 
a5 sles ~ 1684 02774 ° 2429 -0870 eOgi2 
40 ~1824 1549 5220 Beas S) | -0791 SOULE 
45 51670 ellos ~2474 02224 2047 20507 
50 SIS Ome Aes iaiak sees -0341 p04 
5 sLOG2 20850 ~1824 01447 sORT2 AG lesy 
60 NOs: - 0603 02724 R28, - ,.0040 0017 
65 20GL4 - 0206 NOLS 02154 - 0060 S010E9 
Wt. Avg. - 1684 wesc S 5 ORLIES, -LeEG BE TAS) rp bays) 


Table XI. Cross-seasonal comparison of zonally-averaged values of 
Ras a function of latitude. The averaged values for the 2/3-CL 
case are denoted with a sSubScripted 1 and the full-CL case values 
by a subscripted 2. (All values in ly min7t). 
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Mid-Seasonal Dates 


(a) (b) (c) 
16 April 1974 16 July 1974 LGr@etober 1973 
‘ona Sil Ro2 a a7 ui 380 
20S =. 1042 = (USK a esas -.1671 - .0900 -.0942 
LS -.1134 = Ooo -.1772 -.1863 -. 1061 = AO SHS, 
10 =,1297 =.1315 -.1824 =. 1906 -.1107 -. 1093 
> -.1267 -.1284 =-«1L869 ie eZ -.0941 -.0872 
O -.0811 -,1455 -.1604 =, 1663 -.1225 =.1219 
» alloy =,1671 -.1476 =p 527 = 322 -.1334 
10 = 1D7/als} -.1824 -.1522 -, 1582 -.1402 +1434 
dies lS we OS -.1582 -.165/7 -.1488 -, 1495 
20 -.1814 =— os -.1749 -.1853 -.1621 - .1687 
22 -.1823 -.1895 -.1756 -.1846 =i 72 =. 1669 
30 = GOOG -, 1622 =,1673 =e 27) =,17 53 -.1828 
25 —. 1401 esos = 1/757 -. 1824 =, 1613 -.1628 
40 selene 4 =e 253 -.1667 =—en7 Ll -.1707 -.1748 
45 =, 1243 -.1247 =. 1595 -.1612 —Telon -.1718 
50 = 5 JS iLih -.1314 =i 29: -.1277 =, 1639 -.1736 
55 -.1046 = (UII. =O (Ss -.0684 -.1419 -.1453 
60 -.1246 =, JIBS =, L123 -.1046 -, 1486 -.1511 
65 -.0703 = 0 S1E6 -.0648 -.0517 = 1226 =e 2S 
Wt. Avg. -.1457 -. 1486 -.1585 =e leZze -.1450 -.1474 
Table XII. Cross-seasonal comparison of the zonally-averaged values 
of R_- as a function of latitude. The averaged values for the 2/3-CL 


are denoted with a subscripted 1 and the full-CL cases values by a 
subscripted 2. (All values in ly min7}). 


ro 





LAT. 


20S 


138; 


60 


65 


Wt. Avg. 


Table XIII. 


RNMOD 


2/7 3-Cl, 


-.1144 
=eOoZG 
=e Ole is 
-,0554 
-.0042 
-0L0S 
SOE es 
sOo7 2 
Mood? 
~ 9632 
SOB TES 
20958 
~1148 
Oe fey, 
21244 
OCLs 
OZ 


~ 1436 


0288 


Comparison of the ocean-troposphere net radiation R 
k=2, with the model computed net radiation at the top of the 
atmosphere RNMOD and with climatology net radiation at the top of 
the atmosphere from Raschke et al (1973) RN( RAS). Ro and RNMOD are 


Full-CL 


=. 200 
a Lee IL 
-.0844 
sOZes 
-.0345 
(OR 8 
Sy O18 
ae OLOZ 
- 0086 
~0348 
- 0488 
~0546 
jOveZ 
~0504 
~0819 
-Os17 
~ 1034 


- 1082 


BOOZ / 


Shown for both cloud cases. 


Ocean 


Troposphere 

“sil Reo 

= A QUL -,1124 

-,0817 = One, 

20546 =07 72 

- .0468 -.0656 
sO1L3S =O 54 
s02O7 -.0062 
qos ~0145 
20635 cog 
10627 201356 
Jaz wOSc0 
20945 20560 
, LOLS - 0607 
eeZod - 0888 
~0892 Sle A 
~1461 - 1038 
- 1064 - 0766 
- 1603 Bylo leo. 
oO - 1640 
20536 20190 


lll 


RN (RAS) 


=,0504 


me OZ 


OLS 2 


-0616 


~ 1194 


ersa2 


~ 1430 


sLo00 


~ 1426 


- 1484 


~ 1494 


7 Som 


~1164 


0965 


Ooo, 


20715 


20708 


~0858 
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VIILT. ZONALLY-AVERAGED TROPOSPHERIC AND 
OCEANIC HEAT BUDGETS FOR 16 APRIL 1974 
Poe lHE, TROPOSPHERIC HEAT BUDGET 
By averaging over the tropospheric columns for the four meridians 
Pcelayed in Figs. 7, 8, 9 and 10, latitude by latitude, the zonal 
aesteribution of E + Hy results were obtained and are presented in 
Malle XIV. Recall that Roe the tropospheric radiative net cooling is 


listed in Table XII. The resulting heating rate of the troposphere 


May then be expressed as a function of ¢? by the right side of Eq. 8-l 


(ee oe ot (E + Hp) (8-1) 


LAT. 20S 1S 10 5 0 5 10 5 20 
E+H como ome Ole O27 3 720595 6.0592 *<.0609 .1296 -.0/761 
LAT. 25 30 35 40 45 50 oS, 60 65 
E+H aeooee a wwOcor— ole —,049)1 = 70047 —,.0344 —70358 —.0092 —.0303 


Table XIV. Zonally-averaged values of E + Hp as a function of latitude. 
(All values in ly min”). 


S. is the storage heating rate of the column, and oan is the re- 
quired flux-divergence of heat compatible with the heat balance equa- 


tion for the column, for the observed values of Roe Bet Ha sand S.° 


r 
Figure 15 depicts the three heating rates (Ros (E + Hy) and (0 5 43 S_)) 


as a function of latitude for the two cloud-model cases. 
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LATITUDE 


Meoure 15, Latitudinal distribution of R, E + Hype and of oar + S5 
aS a function of latitude for the two 

Ccloud-model cases. Values of E + H, at 15 and 20S were off the 
upper scale of the ordinate and may be found in Table XIV. (Solid 
curves denote 2/3-CL case, dashed curve full-CL). 
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We noted im Section VII.B., the difference Ri - R_»5 is almost negli- 


Grate an this Study, and there 1s little variation as a function of 


PaeLeuldewee tne distribution of E+ He as a function of latitude is also 


Ii 


independent of cloud-model. Therefore, Fig. 15 shows only one curve for 


ae Je The curve of (oe + S ). is shown in Fig. 15, as the residual 


by a’ l 
Siethe right side of Eq. 8-1 for the 2/3-CL case. 


For ease of interpretation in Fig. 15, the curve of Oe + S.)5 Eo 


the full-CL caSe was not plotted. This curve follows closely parallel to 


that Of (Oe + Ss.) being slightly less than (or fo eee thie. latetude 


il al 


range 20S to 45N. From 45N to 65N the curve of (Oe ety S.)5 exceeds the 


curve of (QO. thoes by a small ameunmt. » This crossover effect 1s to be 


al 


expected by examination of the corresponding crossover in the Ri and 


R_>5 curves at these higher latitudes. As is seen in Fig. 15, the net 
heating curve oes oF S. follows the same basic pattern as the E + Hp 
curve, since Ro is relatively constant with latitude. 

The cosine-weighted aii Oe Si -.0512 ly ahne for the average 
tropospheric column. Q., has not been computed in this study; however, 
af Q., + S, is attributed solely to mean storage-cooling S_ of the tropo- 
sphere, the daily storage rate corresponds to a temperature-change rate 
given by 
Pa te = ap 


I’ ° 
Ens x 1440 °C/day (8=2) 


(5) ay 


with AP mb = 800 mb in the troposphere. The resultant cooling rate by 
(8-2) in the troposphere, considering that zero lateral flux-divergence 


applies, is therefore approximately 


eee) et ers hag Os Gh 
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Z 
averaged over the mean cm tropospheric column. 


B. THE LATITUDINALLY-AVERAGED HEAT BUDGET OF THE OCEAN 
The latitudinally-averaged heat-budget of the ocean water-mass 
obeys the following thermodynamic relationship for both the 2/3-CL 


and full-CL cases. 


O pees Ren (Ee A) (8-3) 


vo O r 
Note that Eq. 8-3 is the analog of Eq. 8-1 for the ocean-mass, with 


-(E + H,) representing the average turbulent cooling rate for the 


iA 
water-mass by mixing across the sea-air interface. 

The two terms of the right side of Eq. 8-3 have been computed for 
each meridian and then averaged zonally, with R and E + Hp shown in 
Table XI and XIV, respectively. oS a So has been computed as the 
residual of these two terms at each latitude and is graphed in Fig. 16 
foe both cloud cases. 

The primary ocean-mass heating-function is R, the radiative heating 
rate at the ocean-air interface. As discussed in Section VII, Ry > R, 


for all latitudes considered in this study. Therefore, according to 


ci Oo- 3, SOUS + So) must be greater than cones + S.) by the same 


2 
amounts since -(E + Hp) 1s constant for both cloud cases. This means 
that the internal heating rate of the ocean-mass, aS given in Eq. 8-3, 
1S greater in the 2/3-CL case than in the full-CL case for all latitudes 
considered in this study. On is the required mean oceanic heat-flux 


Givergence at latitude ¢$, and So is the heat-storage in the oceanic 


column. 
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Figure 16. Latitudinal distribution of R, -(E+H.), and of (Oustce) 
as a function of latitude for the two cloud-model cases. The 
value of -(E+H,_) at 15 and 20S were off the lower scale of the 
ordinate and may be inferred in Table XIV. (Solid curves denote 
2/3-CL case; dashed or dot-dashed curves denote full-CL case). 
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The cosine-weighted average of Soi 15 Soy over the latitude range 


of Fig. 16 is 


Sa Sl 
+ =e 
(Q 6 So) 0949 ly min 


(cf£., Tables XI, XIV) which corresponds to a mean heating rate in the 


water-mass column. It should be noted that this heating rate is con- 


Siderably larger (and of opposite sign) than the corresponding atmos- 


pheric effect. 
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IX. CONCLUSIONS 


The present work continues the use of the radiational two-layer 
cloud model adapted from Spaeth's mid-winter study to FNWC gridpoint 
data valid over the four meridians for 16 July 1974. The use of the 
2/3-CL parameterization, as compared to the full-CL after Smagorinsky 
(1960), for the specification of cloud-layer amounts gave the best 
verification for the computed global albedo as a function of latitude 
when compared with observations from NIMBUS III satellite climatology 
comparable to that of 16 July 1974. The use of the 2/3-CL parameteri- 
zation also gave good verification of computed long-wave flux to 
Space when compared to the latitudinal distribution of the same element 
aS given by satellite climatology. 

The chief shortcoming of the computed global albedo versus the 
satellite measured albedo occurred in tropical latitudes (20S - 20N). 
Here, it appears that cloud-elements do not behave as reflecting sheets 
for insolation but rather as focusing devices through the intra-cloud 
spaces. The use of 50-60% model cloud reflectances seems excessive in 
these latitudes. As a practical alternative it may be reasonable to 
tune the cloud reflectances until agreement of the computed global 
albedos with satellite measurements is achieved in the latitude range 
20S - 20N. 

In the topics, all analyses are subject to data-sparsity and hence 
are largely dependent upon climatology and a relatively few radiosonde 


reporting stations. An apparent weakness in this study is that of 
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overspecifying cloud-amounts through the use of Smagorinsky type 
formulas. This could be a result of inadequately resolved relative 
humidities, which in general might well be lower than inferred by 
the FNWC analyses. 

One of the major positive results of the radiative model was 
Eeeeesoimuldtion Of the proper sense of the cross-seasonal radiative 
heat balances at the top of the earth-atmosphere system. A Similar 
positive result applied to cross-seasonal radiative balance at the 


ocean surface centered on 16 July 1974. 
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